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PREFACE 
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Research/ AFSC  under  Contract  F49620-85-C-0013/SB5851-0361  as  a  follow-on  to  work 
done  in  the  summer  of  1987  on  the  Universal  Energy  Systems  Summer  Faculty  Research 
Program  at  Arnold  Engineering  Development  Center,  Arnold  AFB,  Tennessee.  The  author 
returned  to  AEDC  during  the  summer  of  1988  to  continue  work  on  the  project.  Apprecia¬ 
tion  is  extended  to  C.  C.  Limbaugh  and  W.  K.  McGregor  for  suggesting  the  topic,  to  M. 
A.  Simmons  for  running  the  NASA  equilibrium  code,  and  to  Mike  Powell  of  the  Air  Force 
Astronautics  Laboratory  for  providing  numerous  references  to  rocket  test  reports.  M.  K. 
Kingery  supervised  the  Summer  Faculty  program  at  AEDC. 

The  reproducibles  used  in  this  report  were  supplied  by  the  author. 
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1.0  INTRODUCTION 


Many  liquid  rocket  engines  employ  a  film  of  liquid  fuel  as  thermal  protection 
for  the  combustion  chamber  walls.  This  process  was  experimentally  studied  in  the 
1950’s  and  1960's,  however  no  general  analysis  was  developed.  Since  then  a 
number  of  fundamental  heat  transfer  studies  have  been  performed  which  allow  the 
liquid  film  cooling  problem  to  be  analyzed  as  a  collection  of  several  fundamental 
processes. 

The  mechanisms  involved  in  liquid  film  cooling  are  depicted  in  Figure  1. 
Heat  is  transferred  from  the  hot  free  stream  gas  to  the  liquid  film  by  both  radiation 
and  convection.  This  energy  is  absorbed  in  vaporizing  the  liquid  in  the  protective 
film  on  the  wall.  The  vapor  generated  flows  outward  from  the  liquid  film, 
decreasing  the  normally  expected  convective  heat  flux  by  the  well-known 
transpiration  cooling  process.  Downstream  of  the  liquid  film  the  vapor  mixes  with 
the  free  stream  gas  entrained  in  the  boundary  layer,  lowering  the  wall  temperature 
through  the  well-known  "gaseous  film  cooling”  process.  This  provides  thermal 
protection  downstream  of  the  dryout  point. 

The  existing  heat  transfer  correlations  are  for  the  flow  of  a  low  turbulence 
gas  at  constant  velocity  over  a  flat  plate.  These  were  extended  to  the  case  of 
developing,  accelerating,  and  turbulent  free  stream  flows.  Since  the  goal  was  to 
develop  a  simple  model,  suitable  for  design  and  test  analysis,  only  a  simple 
one-dimensional  model  was  considered. 
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2.0  LIQUID  FILM  EVAPORATION 
2.1  CONVECTIVE  HEAT  TRANSFER 

The  flow  over  the  combustion  chamber  walls  is  almost  like  that  over  a  flat 
plate  in  a  uniform  free  stream  flow,  the  difference  being  that  it  arises  from  an 
effective  stagnation  point  rather  than  an  abrupt  leading  edge.  However,  for  a 
turbulent  boundary  layer  the  exact  starting  conditions  quickly  become  unimportant 

In  the  absence  of  the  liquid  film,  the  normal  convective  heat  flux  can  be 
calculated  using  one  of  the  many  flat-plate  data  correlations  (Reference  1).  Most 
agree  within  5%.  One  simple  correlation  with  an  analytical  basis,  termed  "Colburn’s 
Equation",  is  based  upon  a  l/7th  power  law  velocity  profile  (Reference  2). 
Expressed  in  terms  of  the  friction  factor,  Cf,  and  the  Reynold’s  number  based  upon 
the  distance  x  from  the  leading  edge: 

Cfo  =  0.0592  Rex'°  2  (Equation  2.1) 

Using  Reynold’s  momentum-heat  analogy  for  turbulent  convection,  a 
non-dimensional  heat  transfer  coefficient,  the  Stanton  number,  is  expressed: 

Sto  =  v.  Cf0  Pr'°  8  (Equation  2.2) 

Colburn’s  Equation  is  applicable  for  Rex<1.10T.  In  most  rocket  engines  this 
condition  is  satisfied. 

In  a  rocket  engine  there  is  a  large  temperature  difference  between  the 
combustion  gases  and  the  liquid  film,  so  that  the  temperature  at  which  the 
properties  are  evaluated  is  important.  It  is  common  practice  in  boundary  layer 
correlations  to  evaluate  properties  at  the  free  stream  gas  temperature.  A  correction 
factor  for  the  free  stream  to  wall  temperature  ratio  is  then  multiplied,  the  exponent 
being  0.26,  0.4,  or  0.5  depending  upon  the  author  (References  3,  4,  or  5, 
respectively).  However,  evaluating  all  of  the  properties  at  the  mean  temperature  of 
the  free  stream  and  wall  eliminates  such  correction  factors,  and  is  assumed  here. 

As  the  boundary  layer  grows,  the  cylindrical  geometry  becomes  important. 
The  most  extensive  measurements  in  fully-developed  pipe  flow  at  high  temperatures 
were  by  Humble  (Reference  6).  He  obtained  a  correlation  which  can  be  written 
(using  Nu  =  StRePr): 

Sto  =  0.023  ReD'°  • 2  Pr'°  ■ 6  (Equation  2.3) 
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where  the  Reynold’s  number  is  based  on  the  pipe  diameter.  Humble  found  that  if 
the  property  values  were  evaluated  at  the  mean  temperature,  then  no  temperature 
correction  factor  was  necessary. 

Bartz  (Reference  7)  found  that  Equation  2.3  predicted  the  heat  flux  at  the 
throat  of  an  KFNA/N2F4  rocket  if  the  constant  was  0.026  .  In  Bartz’  tests  the  flow 
was  not  completely  developed,  which  explains  the  slightly  larger  constant.  This 
so-called  "Bartz  Equation"  has  been  commonly  used  to  calculate  convective  heat 
fluxes  In  rocket  engines.  Not  surprisingly,  it  usually  underpredicts  the  convective 
heat  flux  upstream  of  the  throat  where  the  assumption  of  fully-developed  flow  is  not 
justified. 


For  x>3-53  D  the  flat  plate  correlation  (Equation  2.2)  predicts  a  lower  heat 
flux  than  that  for  fully-developed  flow  (Equation  2.3).  Certainly  the  heat  flux  cannot 
fall  below  that  for  fully-developed  flow.  To  patch  these  two  limiting  cases  together. 
Equation  2.2  can  be  used,  but  with  x  replaced  by  an  effective  Xg  (Reference  8): 


x 


e 


3.53  D 


_ X 

3.53 


-1/m 


The  exponent  m=  1.2  gives  the  best  comparison  with  the  data  of  Barbin  and  Jones 
(Reference  9),  as  shown  in  Figure  2.  Their  data  was  taken  at  Rej^  =  388,000  which 
is  in  the  range  of  the  data  analyzed  in  Section  4.1.  The  optimum  exponent  is 
probably  a  function  of  the  Reynold’s  number. 


The  convective  heat  flux  for  the  flow  of  free  stream  gas  over  the  dry  walls  of 
the  combustion  chamber  is  then  calculated  as: 


Qconv  =  ^o  ^ 


where  hQ = GfneanCpgStt>  .  Gmean  is  the  gas  mass  flow  per  area,  evaluated  at  the 
mean  temperature  between  die  gas  and  liquid  film.  It  is  scaled  from  the  chamber 
value  (Gch)  as: 

Gmean  =  GchCr8/Tmean)[(Ug-Ul)/Ug] 

The  second  factor  accounts  for  the  velocity  of  the  free  stream  relative  to  the  liquid 
surface.  The  liquid  surface  velocity  (U|)  is  given  in  Section  2.2.2  . 
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2.1.1  Free-stream  Turbulence 

The  convective  heat  flux  increases  by  a  (actor  Kt  =  1  +  4.0  et  in  the  presence 
of  an  rms  turbulence  fraction,  et ,  in  the  free-stream  flow  (Reference  10). 

Two  studies  measured  the  turbulence  intensities  in  rocket  engines.  Hersch 
(Reference  11)  found  et  =  0.10  to  0.05  at  distances  of  2  to  8  inches  from  the 
injector  in  a  LQx/GHa  engine,  while  Talmor  (Reference  12)  found  et  =  0.20  to  0.15 
at  distances  of  6  to  23  Inches  from  the  injector  in  an  NaO<  /Az50  engine. 


2.12  Transpiration  Corrections 

The  total  heat  flux  due  to  both  convection  and  radiation,  Qtot,  is  absorbed  in 
the  liquid  film,  causing  an  initial  temperature  rise:  dT[jq/dx=Qtot/(TCp)  ,  where  r 
is  the  local  liquid  mass  flow  rate  per  circumference.  After  the  liquid  film  reaches 
the  saturation  temperature,  Ty,  the  evaporation  rate  per  area  is: 

•  •  «  • 

®v  =  Qtot/*  =  (Qconv+Qrad) /*  (Equation  2.4) 

This  vapor  flows  away  from  the  liquid  film,  similar  to  gaseous  transpiration 
through  a  porous  wall.  Transpiration  decreases  the  wall  shear  stress  and  convective 
heat  flux  by  a  factor  dependant  upon  the  "blowing  parameter":  F  •  mv/G. 

The  simplest  transpiration  analysis  assumes  that  the  transpiring  vapor  does 
not  change  the  thickness  of  the  viscous  sublayer  in  the  boundary  layer.  Normally 
the  velocity  profile  in  the  viscous  sublayer  is  linear.  With  blowing,  it  becomes  an 
exponential  function.  The  ratio  of  the  wall  shear  stress  in  the  blown  (rw)  to 
unblown  case  (rwo)  is: 

-v  =  — 2 -  where:  z  =  — 

Two  ez  -  1  cfo 

This  analysis  is  termed  the  "Couette  flow  model*  (References  5  and  13)  or 
alternatively  "film  theory*  (Reference  14).  Since  St=Cf/2,  the  convective  heat  flux  is 
decreased  by  the  same  factor.  It  is  convenient  to  rearrange  the  expression  to  a 
logarithmic  form,  where  the  argument  contains  the^  actual  Stanton  number  (St), 
rather  than  the  unblown  value  (StQ).  Substituting  Qconv  =  St  G  Cpg  aT,  gives  a 
simple  result: 

h_  _  lnf.l+H)  where:  H  =  E  =  SpgfaT  +  2radl 

hQ  H  St  X  t  h  J 
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Without  radiation,  the  convective  heat  flux  decrease  depends  only  upon  the 
non-dimensional  factor  CpgAT/i  ,  termed  the  "Spalding  transfer  number"  in  droplet 
evaporation  research.  With  radiation,  h  must  be  determined  implicitly. 

Results  of  more  sophisticated  transpiration  analyses  are  compared  with  die 
Couette  flow  model  in  Figure  3-  Since  the  existing  experimental  data  shows 
considerable  scatter,  there  is  no  advantage  in  using  the  more  complicated  analyses. 

To  account  for  an  injectant  different  than  the  free  stream  gas,  the  parameter 
H  must  be  multiplied  by  a  correction  factor,  Kj^,  being  either  the  vapor  to  gas 
specific  heat  ratio  to  the  power  0.6  (References  5  and  1 5)  or  the  gas  to  vapor 
molecular  weight  ratio  to  the  same  power  (Reference  16),  when  Mc<Mg.  These 
two  forms  are  identical  for  ideal  gases  with  the  same  y.  For  Mc>Mg,  Rubesin 
(Reference  16)  gives  the  exponent  as  0.35.  The  data  of  References  5  and  15  was  for 
Mc<Mg  only  and  thus  does  not  contradict  this  second  exponent  These  corrections 
are  confirmed  by  Landis'  turbulent  boundary  layer  model  (Reference  3). 


Since  the  specific  heats  were  considered  in  the  transpiration  model,  they 
should  not  require  a  separate  accounting.  However,  a  molecular  weight  difference 
would  require  a  correction  term,  because  a  lower  density  injectant  would  displace  a 
greater  volume  in  the  boundary  layer,  decreasing  the  heat  flux.  The  effect  should 
be  less  pronounced  when  the  injectant  has  a  greater  density  than  the  free  stream 
gas,  as  given  by  Rubesin's  second  exponent.  For  this  reason  the  correction  term 
based  upon  the  molecular  weight  ratio,  as  suggested  by  Rubesin,  is  used  in  the 
model  as: 


where: 


a=  0.60 
0.35 


for  Mc<Mg 

for  Mc>Mg 


2.13  Liquid  Film  Disturbances 

Two  types  of  disturbance  to  the  surface  of  the  liquid  film  have  been 
observed.  These  have  been  studied  by  Kinney  (Reference  17),  Knuth 
(Reference  18),  and  Gater  (Reference  19). 

The  first  type  of  disturbance  is  the  appearance  of  random,  small-scale 
structures,  with  a  pebbled  appearance.  These  disturbances  are  always  present  and 
cause  an  effective  surface  roughness  which  may  increase  the  heat  flux. 
Photographic  studies  show  that  the  length  scales  of  these  disturbances  decrease  at 
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higher  velocities,  although  no  general  correlations  were  determined 
(References  17, 18). 

The  other  type  of  disturbance,  realized  at  higher  coolant  flow  rates,  is  the 
appearance  of  large  waves  traveling  in  the  flow  direction.  These  waves  exist  only  at 
the  upstream  positions  where  the  liquid  film  is  thick  enough  to  sustain  them.  The 
mass  loss  rate  in  the  regions  where  these  waves  do  exist  is  2  to  4  times  the  normal 
evaporation  rate,  independent  of  the  film  thickness.  Shearing  of  droplets  from  die 
crests  of  the  waves,  without  evaporation,  is  considered  the  primary  mechanism  of 
increased  mass  loss. 

The  point  of  onset  of  the  large  waves  is  well  defined.  Knuth  was  able  to 
correlate  both  his  data  and  that  of  Kinney  for  the  transition  point.  Assuming  that  the 
mixture  above  the  liquid  him  is  mostly  vapor  and  that  the  vapor  to  liquid  viscosity 
ratio  is  greater  than  0.03,  Knuth’s  correlation  for  the  critical  liquid  mass  flow  per 
circumference  is: 

rcr  =  1.01  x  10B 

Gater  found  quite  different  results.  In  photographic  studies  he  noticed  large 
waves  only  for  water  coolant.  Mass  was  lost  without  evaporation  for  all  of  his  test 
conditions.  Furthermore,  the  mass  loss  rate  was  dependant  upon  the  local  film 
thickness,  proportionate  to  the  local  flow  rate. 

These  differences  might  be  explained  by  the  different  experimental 
arrangements.  Kinney  and  Knuth  both  measured  the  length  of  the  liquid  film  for 
different  liquid  flow  rates  by  determining  the  dryout  point,  using  axially  spaced 
thermocouples.  Gater  used  a  fixed  film  length  and  captured  the  excess  liquid  in  a 
downstream  slot  Obviously  die  potential  existed  for  some  of  the  liquid  to  flow  past 
the  capture  slot.  At  high  injection  rates  it  is  also  difficult  to  insure  that  all  of  the 
liquid  is  placed  on  the  wall.  Knuth  considered  this  problem  in  detail.  Either  of 
these  problems  might  explain  the  different  results  of  Gater.  hi  favorably  comparing 
his  proposed  correlation  with  Kinney’s  data,  Gater  employed  an  adjustable  constant 
and  failed  to  properly  account  for  quantities  on  a  unit  area  basis. 

Liquid  film  waves  have  also  been  studied  with  application  to  chemical 
processes.  Hanratty  and  Hersham  (Reference  20)  give  a  correlation  for  the  free 
stream  velocity  at  the  transition  to  large  waves  in  terms  of  the  gravitational  constant. 
However,  they  were  unable  to  explain  the  results  of  Kinney  and  Knuth  in  a 
horizontal  tube.  Apparendy  their  correlation  is  useful  only  in  vertical  tubes  at  low 
gas  velocities.  Woodmansee  and  Hanratty  (Reference  21)  measured: 
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Tcr=  0-0641  kg/m-s  for  water,  which  is  about  3  times  lower  than  predicted  by 
Knuth’s  correlation.  Of  related  interest,  Tatterson  and  Dallman  (Reference  22)  give 
a  correlation  for  the  mean  diameter  of  the  droplets  sheared  off  of  the  liquid  film. 

The  present  model  is  valid  only  when  the  liquid  flow  rate  is  below  the  critical 
value  given  by  Knuth’s  correlation.  The  small-scale  disturbances  are  not  assumed  to 
increase  the  convective  heat  flux. 


22  RADIATIVE  HEAT  TRANSFER 

In  most  rocket  combustion  chambers  the  radiant  heat  flux  is  negligible  in 
comparison  with  the  convective  heat  flux  (References  23  and  24).  However,  the 
transpiration  of  vapor  from  the  liquid  film  decreases  the  normally  expected 
convective  flux  to  the  extent  that  the  radiant  flux  may  become  dominant  in 
determining  the  liquid  evaporation  rate. 

The  radiation  from  combustion  products  is  difficult  to  calculate  from  first 
principles  due  to  the  complicated  molecular  spectra  of  the  many  species  present. 
Molecular  band  models  attempt  to  replace  the  actual  spectral  lines  present  with  a 
statistical  distribution  of  lines,  having  the  same  gross  properties.  Such  modeling  is  a 
major  research  effort  and  was  not  warranted  in  the  present  study.  However,  with 
access  to  spectral  codes,  it  is  preferred  to  the  simple  approach  described  here. 

Fortunately,  in  a  high  temperature,  high  pressure  environment  many  of  the 
gas  spectral  properties  can  be  more  simply  treated.  High  temperature  causes  each 
spectra]  line  to  be  "doppler  broadened"  due  to  kinetic  motion.  Even  more 
significant  is  "collision  broadening",  due  to  high  pressure.  This  broadening  tends  to 
smear  the  spectrum  lines  together  into  continuous  bands.  'With  a  large  number  of 
different  species,  there  would  be  few  gaps  in  the  spectrum.  Soot  radiation  aids  in 
filling  the  gaps.  In  the  limit,  a  very  dense  gas  mixture  in  thermodynamic 
equilibrium  radiates  as  a  perfect  blackbody. 

22.1  Total  Emit tance  of  Gas 

The  total  emlttance  of  a  gas  is  the  ratio  of  the  radiant  intensity  of  the  gas  to 
that  of  a  blackbody  at  the  same  temperature,  averaged  over  the  entire  spectrum. 
This  quantity  is  much  easier  to  measure  than  the  detailed  spectral  data  and  has  been 
available  for  many  years.  Unfortunately  data  is  available  only  for  H2O  and  CO2, 
however,  this  covers  many  cases  of  combustion  interest  since  symmetric  diatomic 
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molecules,  such  as  Na ,  do  not  radiate  significandy. 


The  most  recent  emittance  data  at  1  atm  pressure  is  compiled  by  Siegel  and 
Howell  (Reference  25).  This  data  is  replotted  in  Figure  4  as  a  function  of  the  optical 
density,  pQpt,  which  is  the  gas  partial  pressure  times  the  path  length  through  the 
gas.  Note  that  at  very  high  optical  densities  the  emittances  reach  limidng  values,  e  f, 
of  0.825  for  H2O  and  0.231  for  CO2.  These  curves  may  be  analytically  fit  to  a 
function: 


*  - e,  [l  +  (*?*]'”  ]"1/n 


with  coefficients: 

TCKI 

H2O 

n 

T  CK) 

CO2 

cfatimml 

n 

1000 

0.165 

0.45 

1000 

0.05 

0.6 

2000 

0.90 

0.  65 

1500 

0.075 

0.6 

3000 

2.05 

0.61 

2000 

0.15 

0.6 

A  three  point  interpolation  is  used  for  temperatures  between  these. 


A  correction  Is  required  for  pressures  other  than  1  atm  by  multiplying  a 
correction  factor,  Kp,  to  each  emittance.  The  curves  for  these  pressure  corrections 
are  fit  by  the  functions: 

H2O:  Kp  =  1  +  Ci (l-exp[ (1-P[1+NW])/CZ] } 

where:  Ci  =  0.26  +  0.74  exp(-2.5  pHao) 

C2  =  0.75  +  0.31  exp(-10  PH20)  atm 

CO2:  log10  Kp  =  0.036  pC02 * 4 * 3  [  1+  ( 2  log10P)-®]-Vm 

where:  n  =  100  pqq2 

An  additional  correction  is  necessary  to  account  for  overlaps  in  the  two 
spectra.  This  correction  (for  T>  1200K)  is  fit  by: 


At  =  0.0551  Kx  [1-exp (-4  popt)]  [1-exp (-12. 5  popt) 1 


where:  Kx  =  1 


2  NW  _  1 
Nw+Nc 


n  =  5.5  [1+ (1. 09  Popt)-8*8®]-1/8**8 
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In  the  equations  above: 


P HsO  »  PC O2  '  Popt  atm.m 
P  <>  atm 

Nw  =  mole  fraction  of  water  in  mixture 
Kc  =  mole  fraction  of  COz  in  mixture 

The  total  gas  emittance  is  then: 

£g  =  £Hso  +  £COa 

To  evaluate  these  terms  it  is  necessary  to  input  the  optical  path  length 
through  the  gas.  Since  this  generally  varies  over  the  sight  angle,  the  most  direct 
method  would  be  to  calculate  the  emittance  at  a  number  of  sight  angles  and  to 
average  them,  weighted  by  the  projected  area  at  each  angle.  However,  a  simpler 
approach  is  to  use  an  overall  effective  length.  A  very  simple  expression  is  quite 
accurate  in  determining  the  effective  length  (Reference  25):  Lgff  =  0.95(4V/A)  , 
where  V  is  the  chamber  volume  and  A  is  the  surrounding  surface  area. 

In  applying  this  formula,  the  downstream  section  of  the  chamber  can  be 
assumed  to  be  an  infinite  cylinder,  for  which  Lgg-  =  0.95  D.  If  the  cloud  of  droplets 
coming  out  of  the  injector  is  assumed  to  be  perfectly  reflective,  then  the  upstream 
direction  can  also  be  considered  an  infinite  cylinder.  Since  the  droplets  have 
numerous  partially  reflecting  surfaces  this  is  probably  the  best  assumption. 

Alternatively,  if  the  droplets  are  assumed  to  be  perfectly  absorbing  then  the 
effective  upstream  length  is:  Leff  =  0.95D[4x/(D  +  4x)].  The  upstream  and 
downstream  emittances  are  then  averaged  to  obtain  the  total  emittance.  With  this 
assumption  the  radiant  heat  flux  must  be  calculated  at  every  axial  position. 

It  is  also  necessary  to  consider  the  possibility  of  reflective  walls.  Reflective 
walls  increase  the  effective  sight  paths.  The  simplest  correction,  due  to  Egleti 
(Reference  26),  is  to  multiply  the  effective  length  above  by  a  factor  Aw  to  the  power 
-0.85  ;  where  Aw  is  the  wall  absorptivity. 

With  the  gas  emittance  determined,  the  radiant  heat  flux  is  calculated  as: 

Qrad  =  °  £  g  ITg<  -  Tv4  ] 

where:  a  =  5-67  10"*  W/ms  .“K* 
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A  final  concern  is  that  the  above  relations  are  for  the  radiation  from  hot  gases 
in  chemical  equilibrium.  However,  there  is  evidence  that  the  radiation  from 
transient  species  can  be  important.  Ziebland  (Reference  26)  measured  an  emittance 
of  0.22  in  a  small  O2  /Hs  rocket  operating  at  10  atm.  The  radiation  peak,  which 
occurred  at  a  position  8  cm  from  the  injector  face,  was  almost  three  times  the 
expected  value.  Some  20  cm  downstream,  the  emittance  settled  to  the  expected 
value.  This  initially  large  radiation  was  attributed  to  radiation  by  transient  OH. 
While  the  radiation  from  OH  in  atmospheric  flames  contributes  only  a  faint  blue 
glow  of  no  significant  intensity,  at  high  pressures  it  becomes  significant 
(Reference  27).  Since  OH  concentrations  are  difficult  to  calculate,  radiation  by  OH 
was  not  considered  in  the  present  model. 


2J.2  Liquid  Film  Burnout 

The  radiant  heat  can  penetrate  the  liquid  film  and  be  absorbed  directly  at  the 
combustion  chamber  walls.  It  is  then  conducted  into  the  liquid  film  by  boiling  heat 
transfer.  With  a  high  enough  heat  flux  the  liquid  film  can  "burnout",  as  in  normal 
pool  boiling. 

The  same  geometry  of  a  thin  liquid  film  flowing  across  a  heated  surface  has 
been  studied  by  Monde  and  Katto  (Reference  28),  Katto  and  Ishii  (Reference  29), 
and  Mudawwar,  et.  al.  (Reference  30).  They  all  correlated  the  burnout  heat  flux  as: 


where  pv  and  p%  are  the  densities  of  the  vapor  and  liquid,  a  is  the  surface  tension, 

X  is  the  latent  heat  of  vaporization,  L  is  the  length  of  the  heated  surface,  and  U  is  the 
average  velocity  of  the  liquid  film. 

The  constants  fit  to  the  data  by  each  group  and  the  parameter  ranges  arc 
given  below: 

Group _ g _ DJ _ na _ U(m/s) _ L(cm) _ 

Monde  and  Katto  0.0591  0.725  0.333  3  to  26  0.56,  1 

Katto  and  Ishii  0.0164  0.867  0.333  1.5  to  15  1-2 

Mudawwar,  et.al.  0.0881  0.867  0.432  0.4  to  2  6.4,  12.7 

Monde  and  Katto  studied  a  circular  jet  of  water  and  Freon®  113  impinging 
downward  on  a  heated  plate.  Katto  and  Ishii  studied  a  plane  jet  of  water, 
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Freon®  11 3,  and  1,1,1-trichloroethane  flowing  across  a  heated  plate.  Mudawwar,  et. 
al.  studied  a  film  of  Flu  or  inert®  liquid  (PC-72)  flowing  downward  along  a  heated 
surface.  Since  they  used  a  single  liquid,  they  assumed  the  same  density  ratio 
exponent,  nj ,  as  Katto  and  Ishii.  Since  Katto  and  Ishli’s  test  conditions  appear  to  be 
closer  to  those  for  liquid  film  cooling  in  rocket  engines,  their  correlation  is  used. 

All  of  the  researchers  observed  that  when  the  critical  heat  flux  was  exceeded, 
the  liquid  film  separated  from  the  heated  surface,  Katto  and  Ishii  observed  that  with 
water,  the  separation  occurred  upon  first  contact  with  the  heated  surface,  whereas 
the  organics  maintained  contact  with  the  surface  for  some  distance  before 
separating.  Mudawwar,  et.  al.  observed  separation  upon  first  contact  with  the 
surface  for  their  Fluorinert®  tests.  They  also  observed  that  the  film  reattached  to 
the  surface  when  the  surface  temperature  was  decreased. 

This  total  separation  of  the  film  from  the  surface  In  the  experiments  may 
have  been  due  to  the  constant  heat  flux  condition  imposed.  In  a  rocket  engine  the 
liquid  film  would  be  less  susceptible  to  burnout  than  in  these  heat  transfer 
experiments  because  separation  of  the  liquid  film  from  the  surface  would  generate 
a  region  of  droplets  and  bubbles  which  would  cause  scattering,  decreasing  the 
radiative  transmission  through  the  film. 

Use  of  these  burnout  correlations  is  questionable  since  they  are  expressed  in 
terms  of  the  overall  heated  length.  Ideally,  the  burnout  point  should  be  expressed 
in  terms  of  local  conditions,  such  as  the  local  film  thickness.  They  are  mainly  useful 
as  an  order  of  magnitude  calculation.  With  this  in  mind,  the  burnout  heat  flux  for 
several  rocket  fuels  is  calculated  below  for  the  case:  L  =  5  cm  and  U  =  1.5  m/s,  at  a 
pressure  of  100  psla  (data  Reference  31): 


Fuel 

Tsat 

(*K) 

X 

(10a J/kg) 

PI 

(kg/ra3) 

Pv 

(kg/m3) 

(10~3N/m) 

Qbo 

(kW/m2 ) 

MMH 

433 

663 

720 

9-79 

20.2 

419 

Az50 

413 

870 

778 

9.50 

18. 

548 

Ha 

29 

335 

70 

5-75 

0.6 

17.6 

In  most  rocket  engines  the  radiant  flux  slightly  exceeds  these  critical  values. 
However,  some  of  the  radiant  heat  is  absorbed  by  the  liquid  film.  Since  hydrogen 
has  such  a  low  predicted  burnout  heat  flux  and  does  not  absorb  Infrared  radiation, 
it  appears  unsuitable  for  liquid  film  cooling. 

Katto  et  al.  suggest  a  correction  factor  when  the  liquid  is  injected  at  a 
temperature  below  the  saturation  temperature,  however  it  was  not  correlated  over  a 
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large  range  and  is  significantly  larger  than  such  correction  factors  for  pool  boiling. 
If  the  burnout  point  is  controlled  only  by  local  conditions,  and  the  liquid  is  heated 
to  the  saturation  temperature  at  die  burnout  point,  then  no  sub-cooling  correction 
should  be  necessary,  therefore  none  is  used  in  the  present  model. 

The  fraction  of  radiation  transmitted  through  the  liquid  film  is:  exp(-ctt)  ; 
where  a  is  the  absorptivity  of  the  liquid,  averaged  over  the  spectrum  of  the 
radiation.  The  liquid  him  thickness,  t,  is  calculated  from  the  wall  shear  stress. 
Assuming  laminar  flow: 

t-p-s-th 

L  p  TwJ 

where  r  is  the  local  coolant  mass  flow  per  circumference.  The  surface  shear  stress, 
tw  Is  calculated  by  Equation  2.1,  with  the  transpiration  correction  of  Section  2.1.2  . 
For  simplicity  h/ho  is  used  in  place  of  tw/two. 

The  average  liquid  film  velocity,  needed  for  the  correlation,  Is:  U  =  T/(pt). 
The  surface  velocity  of  the  liquid,  needed  in  Section  2.1,  is  twice  this  average  Aim 
velocity. 

Using  these  values,  the  burnout  heat  flux  is  calculated  with  Katto  and  Ishii’s 
correlation  at  each  axial  position  in  the  rocket  chamber  and  compared  with  the 
radiant  heat  flux  transmitted  through  the  liquid  film.  A  warning  results  when  the 
radiant  flux  exceeds  the  burnout  heat  flux.  However,  as  discussed  in  Section  4.0, 
this  burnout  warning  was  not  found  to  give  a  reasonable  comparison  with  any  of 
the  data  analyzed. 


23  COMPARISON  WITH  OTHER  ANALYSES 
23.1  Purdue  University  Jet  Propulsion  Center  (JPC): 

A  number  of  liquid  film  cooling  analyses  have  been  presented  by  researchers 
from  the  Jet  Propulsion  Center  at  Purdue.  In  early  studies  Zucrow  and  Graham 
(Reference  32)  attributed  the  increased  mass  loss  rate  after  transition  not  to  large 
surface  waves,  but  to  variations  in  the  termination  point  of  the  liquid  film  around 
the  circumference.  Zucrow  and  Sellers  (Reference  33)  attempted  to  model  the  mass 
transfer  with  no  real  success. 
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In  a  later  analysis  Warner  and  Emmons  (Reference  34)  gave  a  result  which 
was  restated  by  Ziebland  (Reference  26).  For  the  no  blowing  case  it  reduces  to  an 
unusual  result: 

St  =  g  Cf  (PrUg)'1 

It  was  not  possible  to  reproduce  their  calculations  with  the  suggested  constants, 
even  after  correcting  for  an  apparent  sign  error  in  their  Equation  5,  so  their  analysis 
was  not  considered  further. 

Gater,  et.  al.  (Reference  35)  attempted  to  relate  all  of  the  previous 
expressions  from  the  JPC.  In  a  later  report  (Reference  19),  Gater  mentions  the 
transpiration  analysis  of  Section  2.1.2,  but  discounts  it  due  to  the  anomalies  in  his 
data  mentioned  in  Section  2.1.3  ■ 

232  Shembharkar  and  Pai  (Reference  36) 

These  authors  use  a  Prandtl  mixing  length  turbulence  model  for  the  gas 
convection,  coupled  to  a  viscous  flow  model  for  the  liquid  film.  Landis 
(Reference  3)  and  Economis  (Reference  37)  reported  similar  turbulent  boundary 
layer  models  which  gave  no  real  improvement  over  the  simple  Couette  flow 
transpiration  model,  so  the  need  for  this  complexity  is  questionable. 

In  their  datum  case,  the  starting  boundary  layer  thickness  is  12.5  mm, 
corresponding  to  a  downstream  starting  position  of  x=  0.812  m,  by  the  boundary 
layer  growth  relations  of  Section  3-0  .  Since  the  boundary  layer  is  well  developed 
upon  contacting  the  relatively  short  liquid  film,  the  evaporation  rate  should  be 
almost  constant.  Instead,  they  calculate  an  initially  high  evaporation  rate  of 
0.2  kg/s-m2,  dropping  exponentially  to  0.124  kg/s-m2,  possibly  due  to  numerical 
problems.  The  simpler  analysis  of  the  present  report  predicts  an  almost  constant 
evaporation  rate  of  0.106  kg/s-m2 .  The  difference  between  this  and  their  final 
evaporation  rate  is  exactly  accounted  for  by  the  blowing  correction  factor  (h/ho).  It 
is  not  apparent  that  the  momentum  of  the  transpiring  vapor  was  properly  included 
as  a  boundary  condition  in  their  model.  This  vapor  momentum  is  what  causes  the 
convective  heat  flux  decrease  in  transpiration  cooling. 
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2.4  SUMMARY  OF  UQUID  FILM  EVAPORATION  ANALYSIS 

This  section  summarizes  the  equations  and  calculation  order  used  in  the 
liquid  film  evaporation  model.  A  simple  algebraic  formulation  is  given  first,  suitable 
when  radiation  is  negligible  and  the  flow  is  either  fully-developed  or  purely 
boundary  layer. 


24.1  Simple  Expression 

Without  radiation,  the  liquid  film  length,  Lc  ,  is  determined  from  the  coolant 
flow  (per  circumference),  r,  and  the  convective  evaporation  rate,  mconv,  as: 

Lc  =  r/mconv 

Substituting  Equation  2.4,  with  the  transpiration  correction, 

* 

L  =  - E - i - 

c  GCpgAT  StQ(h/ho) 

For  fully  developed  flow  (LC>5D),  St^  is  constant  and  Is  calculated  with 
Equation  2.3,  allowing  1^.  to  be  calculated  above. 

For  pure  boundary  layer  flow  (Lc<  2D),  with  liquid  injection  at  x=  0,  the 
average  Stanton  number  is  1.25  times  higher  than  that  at  x=Lc  (found  by  integrating 
Equation  2.2  from  x=0  tox=Lc).  Thus: 

stolavg  "  1.25[0.0296(GLc/fJg)-°-2Pr“0-6] 


Substituting  above: 

L  =  61-62  u  f _ X _ E - l1*2*  Pro.T6 

C  G  Lcpg  AT  A»g(h/ho)J 


Without  radiation,  the  transpiration  correction  of  Section  2.1.2  reduces  to  a 
simple  form: 

h/ho  =  ln(l+H)  where:  H  =  Ceg^AT  ^ 

H  X 
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2.43  Complete  Formulation 

The  equations  used  in  the  liquid  film  evaporation  FORTRAN  program 
"liquid-for"  are  summarized  below. 

The  radiant  emittance  of  the  hot  combustion  products  is  calculated  first  in 
subroutine  “Emittance*  as: 

Eg  =  Emittance(P,Leff,Tg,Nw,Nc) 

where:  Lgjf  =  0.95  D  Aw‘°  * 

giving:  Qrad  =  o  Aw  eg  [Tg4  -  Tv4  ] 

Subroutine  "Emittance"  is  a  direct  coding  of  the  curve  fits  of  Section  2.2.1. 
Parabolic  interpolation  is  used  between  the  three  temperature  curves  for  C02 
and  HjO  in  function  "Fit".  Since  a  constant  Lgfj  is  assumed  (i.e.  reflective 
droplets),  the  radiative  heat  flux  is  constant,  independent  of  axial  position. 

Axial  steps: 

Beginning  at  the  liquid  film  injection  point,  the  local  convective  coefficient 
without  transpiration  is  calculated  as: 

hQ  =  KtG  CpgSt0 


where:  K*.  =  l+4et  "turbulence  correction  factor" 

StQ  =  ^Cf Pr“° • * 

Cf  =  0.0592  Rex-0 • *  ' 

Rex  =  G  Xe/pg 

G  -  GCh(VTm)[(ug-ul)/ug] 

Xe  -  3 . 53D( 1  +  [x/(3 . 53D) ]—1 •a)° ,l>sa 
Tra  =  \ (Tg+Tv) 

The  convective  heat  flux  with  transpiration  is  calculated  implicidy  as: 

h  =  h0ln(l+H)/H 

where:  H  =  Cpg  KH  /h 

Kh  =  (Mg/Kv)a  i  a=0. 6  if  Mv<Mg 

a=0 . 35  if  Mv>Mg 
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mvap  is  the  liquid  evaporation  rate  per  surface  area,  calculated  below. 

The  wall  shear  stress  is  decreased  by  this  same  factor,  so  that: 

tw  =  two (h/hQ) 

where:  tw  =  *sCfG(Ug-Ui) 

The  convective  heat  flux  is  then:  Qconv  =  h  AT. 

The  Inidal  temperature  rise  of  the  liquid  film  is: 

dTliq/dx  =  Qtot/<rcp)  for  Tliq<Tv 

where:  Qtot~Qconv+Qrad 

Once  Tjjq=Ty,  the  liquid  evaporates  at  a  rate  (per  area): 

•  • 

mvap  "  Qtot/x 

which  decreases  the  liquid  flow  per  circumference  at  a  rate  dr/dx  =  -  m^p. 
The  position  at  which  r  =  0  is  the  "film-cooled  length",  Lg. 

The  liquid  film  thickness  and  average  velocity  are  calculated  using  the  simple 
laminar  "Couette  flow”  result: 

t  =  sqrt[Zpr/(pTw)] 

Uliq  =  r/ot> 

The  surface  velocity  of  the  liquid  is  twice  this  average  velocity:  U|  =  2U|jq 

The  possibility  of  liquid  film  burnout  is  tested.  The  burnout  heat  flux  limit  is 
calculated  as: 

Qbo  =  O.O164X/}<to*S3<pv0'l33(Uliqa/Xliq)0*388 

This  is  compared  with  the  radiant  heat  transmitted  to  the  wall: 

*  • 

Qw  “  Qrad  exp (-at) 

•  • 

When  Qw  >  Qi,0  ,  a  warning  is  printed  in  the  output  file.  However,  this  warning 
has  been  found  unreliable. 
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3.0  GASEOUS  FILM  COOLING 

After  the  liquid  film  has  vaporized  it  continues  to  provide  thermal  protection 
to  the  wall  by  calorimetric  mixing  with  the  hot  free  stream  gas  entrained  in  the 
boundary  layer.  Numerous  correlations  have  been  presented  for  this  process, 
termed  "gaseous  film  cooling".  They  are  normally  expressed  in  terms  of  a  "cooling 
effectiveness": 

r,  -  Tq  ~  Taw 
Tg  -  Tc 

where  Tg  is  the  free  stream  gas  temperature,  Tc  is  the  initial  coolant  temperature, 
and  Taw  is  the  "adiabatic  wall"  temperature  in  the  boundary  layer  which  results 
from  calorimetric  mixing  of  the  two  gas  flows.  The  cooling  effectiveness  is  a 
function  of  a  non-dimensional  distance  (X)  downstream  of  the  injection  point: 

X  =  K  X 

0  4  2  £  •  —1.28 

where:  K  ■  G  Ug  Mc 

• 

Mc  =  coolant  mass  flow  per  circumference 
G  «  Pgug  =  free  stream  mass  flow 
x  =  downstream  distance 


Most  authors  arrange  the  constant  K  in  terms  of  a  coolant  Reynolds  number, 
giving  the  false  impression  that  the  injection  velocity  is  important.  As  long  as  die 
coolant  is  injected  with  low  relative  velocity  so  that  there  are  no  jet  effects,  the  exact 
geometry  of  the  injection  point  is  unimportant. 

The  most  successful  analysis  is  that  of  Kutateladze,  et  al  (Reference  38)  and 
Stollery,  et  al  (Reference  39,  Equation  27)  for  gaseous  film  cooling  on  a  flat  plate  in 
turbulent  flow,  with  coolant  injection  at  the  leading  edge.  To  specify  the  conditions 
of  the  boundary  layer  downstream  of  the  Injection  point,  they  conceptually  replace 
the  injected  coolant  with  an  equivalent  mass  of  free  stream  gas.  To  provide  the 
proper  boundary  layer  growth  rate,  they  identify  an  effective  leading  edge  a  distance 
Xq  upstream  of  the  true  leading  edge  (at  x=0).  Their  analysis  can  be  extended  to 
tiie  case  of  injection  at  a  point  Xj  downstream  of  the  true  leading  edge  as  follows, 
with  the  nomenclature  defined  in  Figure  5. 

From  the  solution  to  the  boundary  layer  integral  equations  for  a  1/7*  power 
law  velocity  profile  (Reference  2),  the  boundary  layer  thickness  and  mass  flow  (per 
width)  a  distance  x*  from  the  effective  leading  edge  is: 
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S  =  0.371x#Rex,-° • 2  where:  x'=  x+x0 

Mjjji  -  7/8  GS  «  0. 325MC[X+X0]  0  *  *  (Equation  3.1) 


The  fictitious  upstream  point  XQ  is  found  from  a  mass  balance  at  the 
injection  point: 


Thus, 


Mbllx=xi  =  Mo  +  »up 

0.325  Mc(Xj+X0)°  *8  =  Me  +  0.325  McX^  ■  *  (Equation  3.2) 
XQ  =  (3.08+Xi° •») 1 *28-  Xi 


where:  =  K  x^ 

The  total  mass  of  free  stream  gas  entrained  up  to  position  x  is: 

•  •  « 

Kg  =  -  Mc  (Equation  3.3) 

From  a  calorimetric  heat  balance: 

fi  +  Cpg  L 
L  CPC  M_ 


7  = 


giving. 


Cpq 

1  =  L1  +  Cpc(°-325lx+xo] 


[“ssiHr 

•*-  »r 


(Equation  3.4) 


Other  correlations  based  upon  similar  analyses  have  been  proposed 
(Reference  39,  40,  41).  However,  most  of  these  are  less  careful  about  accounting 
for  the  effects  of  the  injectant  upon  the  boundary  layer  growth  downstream  of  the 
injection  point.  At  large  downstream  distances  most  reduce  to  the  same  form.  The 
correlation  above  gives  excellent  comparison  with  most  experimental  data 
(Reference  39,  42). 


This  analysis  can  be  extended  to  the  case  of  distributed  injection,  as  exists  in 
the  liquid  film  cooling  process,  with  the  vapor  injected  continuously  along  the 
liquid  film  length.  However  the  resulting  expression  is  unwieldy.  In  any  event,  the 
exact  injection  point  is  of  minor  concern.  Librizzi  and  Cresci’s  analysis 
(Reference  41),  which  gives  predictions  very  close  to  Equation  3.4,  does  not  even 
require  the  position  of  coolant  injection. 


This  analysis  contains  two  assumptions  which  are  known  to  be  wrong.  These 
are  that  all  of  the  gaseous  coolant  remains  in  the  boundary  layer  and  that  all  of  the 
gases  in  the  boundary  layer  are  at  the  same  temperature,  Taw.  In  fact, 
measurements  of  the  concentration  and  temperature  profiles  show  an  S-shaped 
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profile  for  both  (Reference  43).  However,  these  profiles  do  maintain  a  similar 
shape,  growing  away  from  the  wall  at  the  same  rate  as  the  boundary  layer  thickness. 
This  suggests  that  the  above  analysis  can  be  reinterpreted. 

The  S-shaped  temperature  profile  can  be  conceptually  divided  into  a  two 
domain  region,  an  outer  layer  at  the  free  stream  temperature  and  an  inner  layer  at 
Taw.  Since  this  inner  layer  thickness  is  a  constant  fraction  of  the  total  boundary 
layer  thickness,  the  same  mass  ratio  of  hot  gases  will  be  entrained  as  in  the  original 
analysis,  giving  the  same  cooling  effectiveness.  However,  in  this  case  the 
interpretation  is  that  ^aw  is  the  temperature  only  in  the  region  very  close  to  the 
wall. 


Even  with  this  re-interpretation  a  difficulty  remains,  in  that  the  analysis 
predicts  no  free  stream  gas  entrainment  in  fully-developed  flow,  which  is  certainly 
unrealistic  for  turbulent  flow.  Indeed,  Equation  3.4  is  used  in  Section  4.2  for 
fully-developed  flows,  with  satisfactory  results.  That  the  analysis  leading  to 
Equation  3.4  is  apparently  flawed  should  not  detract  from  the  fortuitous  result  that 
it  correlates  all  existing  data  well. 

A  number  of  corrections  to  Equation  3-4  are  necessary  to  account  for  effects 
not  considered  in  the  standard  flat-plate  boundary  layer  analysis.  Corrections  for 
free  stream  turbulence,  foreign  gas  Injection,  thermal  radiation,  and  changing  free 
stream  conditions  are  considered  below. 

3.1  FREE  STREAM  TURBULENCE 

The  flat-plate  boundary  layer  relations  used  in  deriving  Equation  3.4  assume 
no  free  stream  turbulence.  Two  experimental  studies  measured  the  effects  of 
turbulence  upon  gaseous  film  cooling  (References  44  and  45).  A  possible 
correction  to  the  above  analysis  is  to  multiply  the  non-dimensional  distance  X  by  a 
factor  Kt  -  l  +  Cj.  This  allows  the  increased  mixing  caused  by  the  free  stream 
turbulence  to  be  accounted  for  as  an  effectively  larger  downstream  distance.  When 
applied  to  the  data  of  Marek  and  Tacina  (Reference  44),  Cj  is  found  to  be  constant 
with  X  and  varies  with  the  free  stream  turbulence  rms  fraction,  as:  Cj  =  8.67  et , 
for  et  =  0.07,  0.14,  and  0.23  -  For  et  =  0.35  ,  the  data  follows  this  relation  for 
X<  1.5,  thereafter  Cj  increases  as  X°  •* . 

Unlike  Marek  and  Tacina,  Carlson  and  Talmor  (Reference  45)  did  not 
direedy  measure  die  turbulence  level,  but  inferred  it  based  on  screens  placed  in  the 
flow.  The  length  scales  of  their  turbulence  may  also  have  differed  considerably 
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from  Marek  and  Tacina's.  Analyzing  their  data  gives:  Cj  =  1 1-7  et .  In  their  data  Cj 
varied  slightly  with  X  at  all  turbulence  levels.  The  expression  above  is  an  average 
over  X.  Since  their  greatest  turbulence  level  was  22%,  the  increased  mixing  found 
by  Marek  and  Tacina  at  35%  turbulence  was  not  confirmed.  For  the  present  model 
die  two  results  were  averaged,  giving  a  correction  term: 

Kt  -  1  +  10.2  et 

32  FOREIGN  GAS  INJECTION 

Goldstein,  et.  al.  (Reference  43)  found  that  the  effectiveness  values  were 
about  30%  higher  than  expected  when  helium  replaced  air  as  the  coolant  injected 
into  a  heated  airstream.  These  results  were  confirmed  by  Burns  and  Stollery 
(Reference  46).  However,  Carlson  and  Talmor  (Reference  45)  questioned  whether 
the  effect  should  be  attributed  to  the  coolant  properties,  since  the  helium  injection 
velocities  were  lower  than  the  air  injection  velocities  at  the  same  effectiveness 
values. 


In  the  present  model  the  effect  is  assumed  to  be  real  and  is  attributed  to  the 
molecular  weight  difference,  for  the  same  reasons  stated  in  Section  2.1.2  .  Goldstein 
fit  his  data  with  an  empirical  equation  (Reference  43,  Equation  12).  Apparently  he 
failed  to  recognize  that  this  empirical  equation  could  be  put  into  the  form  of 
Equation  3.4  by  multiplying  a  correction  factor,  Kjyj,  of  0.76  to  the  (Cpg/Cpc)  ratio. 
Assuming  a  power  function  dependance  of  this  correction  factor  upon  the 
molecular  weight  ratio,  and  generalizing  from  this  specific  case  of  a  coolant  to  gas 
molecular  weight  ratio  of  0.138,  the  correction  factor  can  be  written  in  general  as: 

Km  =  (Mc/Mg) 0  * 1  * 

With  these  corrections  for  turbulence  and  foreign  gas  injection,  Equation  3-4 
becomes: 

1  -  [l  *  KM(g|}[0'325(KtX+Xo>°-‘  -  ‘j]"1  (Equation  3.5) 
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33  THERMAL  RADIATION 

The  radiant  heat  flux  is  transmitted  through  the  boundary  layer  and 
absorbed  at  the  wall.  This  heat  is  then  conducted  back  into  die  boundary  layer 
gases,  requiring  a  wall  temperature,  Tw»  in  excess  of  the  boundary  layer 
temperature,  Taw: 

Tw  =  Taw  +  Qrad/h 

This  form  is  suggested  for  gaseous  film  cooling  with  a  non-adiabatic  wall 
(Reference  47)  The  heat  transfer  coefficient  Is  calculated  by  Equation  2.2, 
evaluated  at  x’  =  x+xo  . 

The  radiant  heat  also  enters  into  the  heat  balance  of  the  boundary  layer 
gases.  Assuming  that  it  is  spread  evenly  over  the  entire  boundary  layer  mass,  the 
rate  of  temperature  increase  with  distance  is: 

• 

ATaW  =  Qgad 
AX  MblCpg 


3.4  NON-UNIFORM  FREE  STREAM  FLOW 

Equation  3.5  was  derived  for  a  constant  free  stream  gas  flow  rate  and 
temperature.  If  these  quantities  change  as  the  flow  proceeds  downstream,  then 
Equation  3-5  cannot  be  directly  applied.  To  demonstrate  this,  consider  a  sudden 
decrease  in  the  flow  rate  per  area,  G.  Applying  Equation  3-5  locally  would  imply  a 
sudden  decrease  in  the  non-dimensional  distance,  X,  and  a  step  increase  in  the 
effectiveness,  i? .  This  would  suggest  that  the  gases  spontaneously  unmix,  violating 
the  second  law  of  thermodynamics. 

Due  to  the  parabolic  nature  of  the  boundary  layer  equations,  the  local 
conditions  determine  only  the  rate  of  change,  in  the  downstream  direction,  of  the 
boundary  layer  properties.  Equation  3.5  is  an  integral  result  for  the  special  case  of 
constant  flow  conditions.  However,  it  can  easily  be  rearranged  in  the  form  of  a 
differentia]  equation. 

The  rate  of  change  of  *p  with  position  x  is  determined  by  differentiating 
Equation  3-5  .  All  parameters  are  held  constant  in  differentiating,  since  such  was 
assumed  in  deriving  the  integral  expression.  Expressing  the  result  in  terms  of  the 
local  value  of  ij : 
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-  -0.1963  K  «*[«■]  t*[l  -  JgjfJ  -1]]'' 

tj  is  found  by  numerically  integrating  this  egression.  The  initial  condition 
on  ij  at  the  injection  point  could  be  found  from  Equation  3-5.  However,  this 
predicts  ij  <1  at  the  injection  point  when  injection  is  downstream  of  the  leading 
edge,  which  is  unsupported.  A  more  realistic  assumption  that  jj  =  1  at  the  injection 
point  is  used  as  the  initial  condition  in  the  present  study.  The  above  approach  is 
termed  the  "differential  t?"  formulation  to  distinguish  it  from  the  "integral 
correlation"  of  Equation  3.5  . 

The  free  stream  mass  flow  rate  per  area,  G,  varies  inversely  with  the  local 
cross-sectional  area  as: 

G  =  GchCAch/A) 

Increases  in  G,  as  the  nozzle  converges,  increases  die  rate  of  free  stream  gas 
entrainment  through  the  coefficient  K. 

The  other  effect  of  free  stream  acceleration  is  a  drop  in  the  static  temperature 
of  the  free  stream  gas.  Upon  entrainment  into  the  boundary  layer  the  original 
stagnation  temperature  is  recovered  as  the  kinetic  energy  is  converted  back  into 
thermal  energy.  However,  some  of  this  heat  is  then  lost  by  conduction  back  into  the 
cooler  free  stream  gas.  The  result  is  that  the  free  stream  gas  attains  a  temperature 
in  the  boundary  layer  lower  than  the  stagnation  temperature.  This  is  termed  the 
"recovery  temperature",  Tr,  and  is  normally  a  constant  fraction  of  the  difference 
between  the  static  and  stagnation  temperatures: 

Tr  =  T0-(l-r) (T0-Ts) 

where  T0  and  Ts  are  the  static  and  stagnation  temperatures  of  free  stream  gas  and  r 
is  the  "recovery  factor". 

This  recovery  temperature  should  replace  Tg  in  the  equations  above.  For  a 
flat-plate  boundary  layer  without  gaseous  film  cooling,  r  =  Pr1  /3  (Reference  48).  It 
is  assumed  that  this  relation  holds  for  the  gaseous  film  cooling  case,  although  this 
assumption  is  questionable. 

From  the  isentropic  relations  for  compressible  flow: 

Ts  =  T0/[l+  K(7-1)M*] 
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The  local  Mach  number,  M,  is  found  implicitly  from  one 

M  =  A  r_2_  +  Mipr/2 
At  l_7+i  7r] 

M  =  |\r(MAt)2/7r  -  55 


of  two  forms: 


for 

M<1 

for 

M>1 

where : 


At  -  (A/Athroat) 


Tr 


=  adbl 

7-1 


3.4.1  Differential  Entrainment  Formulation 

When  the  recovery  temperature  changes,  the  "differential  17"  formulation 
becomes  invalid,  for  the  subtle  reason  that  the  local  value  of  17  no  longer 
determines  the  local  state  of  the  boundary  layer.  This  problem  is  remedied  by 
recognizing  that  the  true  local  state  of  die  boundary  layer  is  determined  by  the  local 
mass  flow  rate  in  the  boundary  layer,  My.  From  Equation  3.1,  the  local 
entrainment  rate  is: 

-  0.1963  G  |j|g-J"0*a5  (Equation  3.6) 

The  calorimetric  heat  balance  in  differential  form  is: 

d^bl  cpg  ( Tr  ~Taw )  =  1  Mgcpg+McGpc )  dTaw 

The  specific  heat  on  the  left  side  should  be  evaluated  at  the  mean  temperature 
of  Tr  and  Taw,  while  those  on  the  right  should  be  evaluated  at  Taw.  Substituting 
Equation  3-3,  and  dividing  both  sides  by  dx: 

dTaw  =  dMfri  (Tr-Taw)  f^bl  +  f Cpc  -  l]Mcl-^"  (Equation  3.7) 

dx  dx  e  L  lcP9  J  J 

To  apply  these  relations,  Equations  3-6  and  3-7  are  integrated  to  determine 
My  and  "^aw  at  each  x  position.  As  stated,  Taw  is  equated  to  Tc  (or  f?  =  1)  at  x=0. 
The  initial  condition  on  My  is  found  from  Equation  3-2: 

Mbilxi  “  Hc[l+0.325Xi0-8] 

In  applying  this  to  the  curved  walls  of  a  rocket  chamber,  x  is  measured  along  the 
contour  of  the  nozzle. 
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This  "differential  entrainment"  formulation  is  the  most  general.  Although  the 
derivation  is  more  involved,  it  requires  no  more  computation  than  the  previous 
"differential  q”  formulation.  For  the  case  of  constant  free  stream  gas  temperature, 
where  the  local  value  of  q  is  directly  related  to  both  formulations  are 
equivalent.  And,  of  course,  for  constant  free  stream  conditions,  both  reduce  to  the 
integral  correlation  (Equation  3 .5). 

To  include  the  effects  of  free  stream  turbulence  and  foreign  gas  injection,  the 
coefficient  in  Equation  3.6  is  multiplied  by  the  factor  Kf  ,  and  the  (Cpc/Cpg)  ratio  in 
Equation  3-7  is  multiplied  by  1/Kj^.  Both  corrections  are  consistent  with  the 
placement  of  these  factors  in  the  differentia]  q  formulation. 

3.4.2  Circumferential  Change: 

In  a  rocket  engine  the  contraction  of  the  nozzle  decreases  the  circumference, 
causing  an  increase  in  the  boundary  layer  mass  flow  per  circumference: 

dMbl  =  _  Mbl  1  dD 
dx  c  D  dx 

Since  this  increase  is  not  due  to  free  stream  entrainment,  it  does  not  enter 
into  the  energy  balance.  However,  it  is  important  in  that  it  affects  the  total  boundary 
layer  flow  rate,  which  is  found  by  adding  the  changes  due  to  free  stream 
entrainment  and  this  circumferential  change.  This  total  boundary  layer  flow  rate 
determines  the  local  rate  of  free  stream  gas  entrainment,  as  given  by  Equation  3.6  . 
In  adding  this  term  to  the  analysis,  the  coolant  flow  rate  per  circumference  should 
also  be  calculated  locally,  by  scaling  the  chamber  value. 


3.5  TURNING  EFFECTS 

Carlson  and  Talmor  (Reference  45)  measured  the  effectiveness  downstream 
of  a  sharp  turn  in  a  rectangular  duct.  Their  results  can  be  fit  as: 


q  =  [  1  +  CX°  )'« 

with  coefficients: 

turning  angle  C _ & 


30  * 

0.40 

2.2 

45 " 

0.75 

2.1 

60' 

1.40 

2.0 
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To  reduce  approximately  to  Equation  3.4  the  exponent,  n,  should  approach 
0.8  as  the  turning  angle  approaches  0° .  Since  the  exponent  is  increasing  in  the 
opposite  direction  their  results  are  somewhat  anomalous.  For  unturned  flows  in  a 
circular  duct  their  results  followed  Equation  3-4.  Possibly  these  results  were  due  to 
disturbance  of  the  boundary  layer  by  the  sharp  turn  in  their  wind  tunnel. 

Ewen  and  Rousar  (References  49-51)  studied  gaseous  film  cooling  in 
converging  nozzles.  Their  data  is  analyzed  in  Section  4.0  and  leads  to  a  correlation 
for  an  increased  entrainment  in  terms  of  a  "centrifugal  parameter". 


3.6  CHAMBER  CONTOUR 


In  general  the  shape  of  most  rocket  combustion  chambers  can  be  expressed 
in  terms  of  a  few  parameters:  1, ,  rt ,  c,  Dt,  r2 ,  and  0d,  as  shown  in  Figure  6.  Given 
these,  the  limits  of  each  segment  are  expressed: 


la  =  lt  +  Tj  sin0c 

Da  =  Dch  "  2  ri (l-cos*c) 

D*  =  °t  +  2  *2(l-cos0c) 

ls  =  1*  +  (D2-D3)/(2tan*c) 
lt  ■  1,  +  r2  sin0c 
1*  =  It  +  ra  sinfld 
Db  =  Dt  +  2  r2(l-cos0d) 

The  chamber  diameter  is  then  calculated  at  any  axial  position,  1,  by: 


For:  lj  <  1  <  12 

12  <  1  <  13 

ls  <  1  <  16 


D  =  Dch  -  2(r1-sqrttrl2-(l-l1)2]) 
D  =  Da  -  2 (l-la)tanflc 
D  =  Dt  +  2 (r2-sqrt [r2  2- (l-ls ) 2 ] ) 

D  =  Ds  +  2(l-lB)tan0d 


as: 


Is  <  1 

The  downstream  position  along  die  wall,  x,  is  related  to  the  axial  position,  1, 
(fix)2  =  (fil)2  +  (fiD/2)2 


This  expression  is  integrated  as  axial  steps  fil  are  made,  to  determine  the  contour 
length. 
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3.7  SUMMARY  OF  GASEOUS  FILM  COOLING  ANALYSIS 

This  section  summarizes  the  equations  and  calculation  order  used  in  the 
gaseous  film  cooling  FORTRAN  computer  program  "ga5.fbr".  A  simple  formulation 
is  given  first,  suitable  when  the  free-stream  mass  flow  and  temperature  are  constant 
and  radiation  is  negligible. 


Simple  Expression: 


With  constant  free  stream  conditions,  the  wall  temperature  is  easily 
calculated  from  Equation  3-5: 

Tav  -  Tr  “  n  (Tr-Tc) 
where : 

1}  -  |^1  +  Kjjj^ESj  |^0.325 {KtX+Xo}°  *•  -  1  j  j 


-1 


X  =  K  X 

Q  .35  •  -1-25 

K  -  G  nq  Mc 


Complete  Formulation: 


Starting  at  the  end  of  the  liquid  film,  steps  dx  along  the  wall  contour  are  made.  The 
increase  in  boundary  layer  flow  due  to  free  stream  entrainment  is: 

dM,  =  0.1963KtG(pg/Mbl|i)-°«25  dx 

where  the  turbulence  correction  Kf  =  1  +  10.2et  .  The  new  boundary  layer  flow  is 
projected  as:  Mblli  +  1  =  My  |  j+dMt .  Following  the  2nd  order  Runge-Kutta 
method,  the  change  is  recalculated  (dM3)  using  My  |  j+  j  .  The  average  of  the  two 
is  used: 

•  •  • 

dMe  =  h (dMt  +dM2 ) 

The  increase  in  boundary  layer  flow  (per  circumference)  due  to  chamber 
contraction  is: 

dMc  “  -  Mbl{l/D) (dD/dx)dx 

where  dD/dx  is  the  rate  of  change  of  diameter  with  contour  length  x,  as  expressed 
in  Section  3-6. 
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The  boundary  layer  flow  at  the  next  axial  step  is  then: 

Mbl!i+1  “  Mbl| i+(dMe+dMc) 

The  increase  in  boundary  layer  temperature  due  to  the  free  stream  mass 
entrainment  is: 

dTe  =  dMe  (Tr-Taw)[Mbl  +  Mc  (1m  CBS  _  ijj-1 

The  increase  in  boundary  layer  temperature  due  to  thermal  radiation  is: 

d^rad  =  t^rad/ CcpgMbl) 

The  boundary  layer  temperature  at  the  next  step  is  then: 

dTawli+l  "  dTawli+<dTe+dTrad) 

The  wall  temperature  is  higher  than  the  boundary  layer  gas  as: 

Tw  =  Taw+Qra{j/h 

where  h  is  the  convective  heat  transfer  coefficient,  as  determined  in  Section  2. 1  . 

The  recovery  temperature  (Tr)  is  determined  by  the  local  area  ratio,  as 
discussed  In  Section  3  4  .  The  local  area  ratio  is  determined  by  the  contour 
equations  of  Secdon  3.6  . 
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4.0  COMPARISON  WITH  EXPERIMENTS 

4.1  LIQUID  FfLM  EVAPORATION 

4.1.1  G.R.  Kinney  et  al  air-water  experiment  (Reference  17) 

These  heat  transfer  measurements  were  performed  at  fairly  low  temperatures 
in  a  tube  with  almost  folly-developed  flow.  This  made  radiation  negligible  and  the 
evaporation  rates  almost  constant  over  the  liquid  film  length.  For  these  reasons  the 
evaporation  rate  could  be  calculated  by  hand  using  the  formulation  of  Section  2.4.1. 
For  completeness  the  computer  program  was  used  to  include  the  slight  effects  of 
radiation,  interface  motion,  and  boundary  layer  development 

The  static  pressure  varied  from  1.4  to  2.47  atm  over  the  runs.  An  average 
saturation  temperature  of  240  F  is  stated  by  Kinney,  corresponding  to  a  saturation 
pressure  of  1.7  atm.  However,  a  liquid  interface  temperature  of  200  F  is  given  in 
Kinney's  Figure  4.  This  difference  might  be  attributed  to  mass  transfer  to  the 
unsaturated  gas.  A  pressure  of  1.7  atm  and  an  interface  temperature  of  200  F  was 
assumed  in  all  of  the  calculations,  since  no  itemization  was  given  for  each  run. 

The  free  stream  gas  consisted  of  the  products  of  gasoline/air  combustion. 
For  property  values,  pure  air  was  assumed.  For  radiation  calculations  10%  C02 
and  10%  HaO  were  assumed  for  all  runs.  In  any  event,  radiation  was  insignificant 
at  die  low  temperatures  of  the  tests.  Since  the  free  stream  turbulence  was 
unknown,  a  value  of  0%  was  used.  The  property  values  common  to  all  runs  were: 

Gas:  Coolant: 

air  water 

M=  29.  M=  18.  /i=  3.03  I0“4kg/m-s 

P=  1.7  atm  Tv-  366  K  962  kg/m3 

Tc=  300  K  pv=  1.01  kg/m3 

Cpl=  4210  J/kg-K  a=  60.3  10"3N/m 
X=  2.27  10® J/kg  a=  0. 
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The  property  values  of  air  are  listed  below  (Reference  54)  for  each  free 
stream  gas  temperature  used  in  Kinney’s  tests. 


Tg 

(F) 

800 

900 

1200 

1400 

1600 

Tg 

(K) 

700 

756 

922 

1033 

1144 

Pq 

(kg/m3) 

0.857 

0.794 

0.651 

0.582 

0.525 

*  g 

(10“skg/ia-s) 

3.349 

3.515 

3.976 

4.264 

4.550 

at 

Trnean  W  “ 

533 

561 

644 

700 

755 

Cpg 

(J/kg-K) 

1036 

1042 

1062 

1075 

1088 

(10-skg/m-s) 

2.80 

2.90 

3.17 

3.35 

3.51 

Pr 

0.698 

0.698 

0.701 

0.702 

0.703 

Kinney  stated  only  the  free  stream  Reynold’s  numbers.  Using  the  viscosity 
values  above,  the  gas  mass  flux  (Gch)  for  each  test  was  back-calculated  from  the 
stated  Reynold’s  number.  The  liquid  film  length  (1^)  at  an  average  coolant  flow 
r=0.08kg/s-m  was  determined  for  each  test  and  the  average  evaporation  rate 
calculated  as:  mv  =  T  /L^, 


Similarly,  the  experimental  evaporation  rates  for  both  normal  evaporation 
(mVgxp  1 )  and  after  transition  to  the  large  waves  (mvexp2)  were  calculated  from  the 
slopes  of  the  data  plots  in  Kinney’s  Figure  6,  giving: 


TEST 

ReD 

Gch 

LC 

mvcalc 

mvexpi 

mvexpj 

(105) 

(kg/s-m2) 

(cm) 

l 

kg/s-m2 

3 

2  in  smooth: 

800F 

4.4 

290.1 

86.6 

0.0924 

0.141 

— 

800F 

5.6 

369.2 

71.3 

0.112 

0.200 

— 

1200F 

3.5 

273.9 

50.8 

0.157 

0.236 

0.450 

1600F 

2.5 

223.9 

38.5 

0.208 

— 

0.611 

4  in  smooth: 

9  OOF 

6.0 

207.6 

106.7 

0.0750 

0.0794 

0.326 

900F 

8.2 

283.7 

83.1 

0.0963 

0.109 

0.226 

900F 

9.9 

342.5 

71.5 

0.112 

0.153 

0.366 

1400F 

5.3 

222.4 

54.7 

0.146 

— 

0.440 

1600F 

4.7 

210.5 

48.0 

0.167 

— 

0.629 

4  in  rough: 

800F 

6.8 

224.1 

119.9 

0.0668 

0.131 

— 

800F 

11.0 

362.6 

81.5 

0.0982 

0.226 

0.458 

1400F 

5.4 

226.6 

53.9 

0.149 

0.326 

0.930 

1400F 

6.7 

281.2 

45.4 

0.176 

— 

0.896 
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Comparison  between  the  calculated  and  experimental  (mvexp  t )  values  is 
given  in  Figure  7.  While  die  calculations  do  not  compare  favorably  with  all  of  the 
measurements  in  an  absolute  sense,  they  do  correlate  the  data  well.  The  calculated 
evaporation  rates  are  a  factor  35%  lower  than  the  data  for  the  two  "smooth"  tubes. 
Assuming  a  free  stream  turbulence  Intensity  of  8.75%  would  account  for  this 
discrepancy,  as  discussed  in  Section  2.1.1. 

Several  other  factors  could  contribute  to  the  discrepancy.  Most  significant  is 
that  the  actual  propellant  flow  rates  were  not  stated,  so  that  the  gas  mass  flux  (Gcj,) 
had  to  be  back-calculated  without  knowing  the  value  of  gas  viscosity  Kinney  used  in 
calculating  the  stated  Reynold's  numbers.  Also,  the  specific  heat  of  the  actual 
combustion  products  is  slightly  higher  than  that  for  the  pure  air  assumed.  Finally, 
the  "pebbling”  of  the  liquid  surface  observed  by  Kinney  would  give  an  effect  similar 
to  pipe  roughness,  increasing  the  convective  heat  flux. 

The  model  is  not  able  to  predict  the  evaporation  rates  when  large  waves  are 
present,  given  by  the  mvexp2  values.  When  the  large  waves  are  present,  the  rates  of 
mass  loss  are  increased  by  a  factor  of  2.56,  on  average. 
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4.1.2  E.L.  Knuth  air-water  tests  (Reference  18) 


These  tests  were  performed  shortly  after  Kinney  in  a  2.9  inch  diameter  tube 
with  fully-developed  flow.  The  gas  properties  are  again  taken  as  those  of  pure  air. 
10%  C02  and  10%  H20  in  the  products  is  assumed  for  the  radiation  calculations. 
The  liquid  interface  temperature  is  taken  as  339  K,  based  upon  Knuth's  Figure  11. 
This  disagrees  with  the  stated  static  pressure  of  1  atm.  The  discrepancy  is 
attributed  to  mass  transfer  to  the  unsaturated  gas.  The  common  properties  used 
were: 


Gas; 
air 
M=  29. 

P=  1  atm 


Coolant; 

water 

M=  18.  4.26  10“4  kg/rn-s 

Tv=  339  K  P!=  980  kg/ra3 

Tc=  300  K  pv=  0.1664  kg/m3 

Cpl=  4188  J/kg-K  o=  65.1  10~3N/m 
1=  2.34  10* J/kg  a=  0. 


The  gas  properties  used  in  each  calculation  were  (Reference  54): 


20-30: 

45-54 ; 

89-99 

Tg  (K) 

613 

901 

1230 

P g  (kg/m 3) 

0.576 

0.392 

0.287 

at-  Tmean  (K)  — 

476 

620 

785 

Cpg  (J/kg-K) 

1025 

1056 

1095 

fiQ  (E-5kg/m-s) 

2.59 

3.10 

3.60 

Pr 

0.699 

0.700 

0.703 

With  these  property  values,  the  liquid  film  lengths  and  evaporation  rates  for 
each  test  were  determined  at  a  coolant  flow  of  0.08  kg/s-m.  Unlike  Kinney,  Knuth 
gave  the  mass  flow  of  the  freestream  gas  directly,  eliminating  one  unknown. 


•  • 


TEST 

G 

(kq/s-m2) 

Lc 

(cm) 

mvcalc 

(kq/s-m2) 

mvexp 
f kq/s-m 

20-30 

181.0 

165.9 

0.0482 

0.0586 

45-54 

91.6 

125.3 

0.0638 

0.0586 

89-99 

69.2 

92.3 

0.0867 

0.0633 

The  calculated  values  agree  fairly  well  with  the  measured  values,  however 
there  is  not  enough  data  to  generalize. 
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4.13  G.  Morrell  LOx/NHa  rocket  (Reference  52) 

These  test  were  performed  in  a  4  inch  diameter  rocket,  with  liquid  oxygen 
and  ammonia  propellants.  The  liquid  coolant  was  injected  2.8  inches  downstream 
of  the  propellant  injector.  Water,  ethanol  and  ammonia  were  tested  as  liquid 
coolants.  The  ammonia  was  at  super-critical  conditions  and  no  sharp  evaporation 
points  were  determined,  hence  it  is  omitted  here.  Since  the  free  stream  turbulence 
was  not  measured,  a  value  of  zero  was  used  in  the  calculations. 

The  pressure  varied  from  16.9  to  17.7  atm  in  the  water  tests.  For  simplicity 
the  following  saturated  coolant  properties  at  17.4  atm  were  used  in  all  calculations: 

coolant  Tv  X  Cpi  pv  °  Pi  Pv  rcr 

(K>  fio,fel  [»Ll  fiB-sl  f10  ‘SI  f-5  1  fcl 

water  480  1.91  4530  15.9  129  36.2  857  9.09  0.198 

ethanol  448  0.574  5126  13.1  121  7.91  637  27.6  0.143 

The  properties  for  water  are  from  Reference  54.  Most  of  the  ethanol  properties  are 
from  Reid  (Reference  55):  Tv  is  from  Appendix  A,  X  from  Equation  7-9.4,  p\  from 
Equation  3-11.9,  fiv  from  Table  9-2,  and  /q  from  Table  9-8.  The  remaining  ethanol 
properties  (Cpj,  a,  and  pv)  are  from  various  tables  in  Perry’s  Handbook 
(Reference  56). 

The  free  stream  gas  differed  in  composition  for  each  run,  with  the  O/F  mass 
ratio  given  by  Morrell.  The  water  coolant  tests  were  all  lean,  with  a  reaction: 
x,NHa  +  x2Oa  -  1.5xtH20  +  +  (x3-0.75xt)O2 

The  water  vapor  fraction,  xHjq,  was  related  to  the  O/F  mass  ratio  by: 

xH20  =  48/ [40+17 (O/F)] 

The  molecular  weight  of  the  products  was  calculated  as: 

Mg  _  544  f 1+  fO/F) 1 
40+17 (O/F) 

The  following  properties  were  used  for  each  product  species  in  the  gas,  at  a 
mean  temperature  of  1720  K: 


HjO: 

N  2  : 

o2  : 


CP 

(J/k5— K)„ 


Vq 

(10~*  Rcr/m-sl 


Pr 


2684 

1266 

1167 


5.85 

5.61 

6.78 


0.9 

0.708 

0.75 
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The  specific  heats  were  calculated  from  Table  3-2  of  Holman  (Reference  57) 
for  low  pressure.  From  Holman’s  Figure  7.8,  the  corrections  for  chamber  pressure 
are  negligible.  The  viscosities  were  calculated  from  the  Chapman-Enskog 
correlation,  using  die  values  in  Appendix  B  of  Reid  (Reference  55).  The  Prandd 
numbers  were  extrapolated  from  the  charts  in  Appendix  A  of  Kays  (Reference  5). 

The  viscosity  of  the  products  was  calculated  at  an  O/F  mixture  of  1.61,  using 
the  mixture  weighting  of  Wilke  (Reference  55,  Eq.  9.5-13),  giving:  fig  ~  5.8 6  10*5 
kg/m-s.  The  specific  heats  were  weighted  by  the  mass  fractions  of  the  product 
species  and  the  Prandd  numbers  by  the  mole  fractions  for  each  run.  With  these 
property  values  the  liquid  film  lengths  were  calculated  for  each  test: 


Test 

O/F 

T9 

m 

p 

fatml 

G 

fJsg_lf 

1  s-m2  f  1 

cpg 

J  1 

ker-Kl 

Pr 

1 

r 

\te_\ 

s-mi 

Llmeas 

fcml 

Ll calc 

_ fcml 

water 

1 

*  * 

» 

1.61 

2950 

17.4 

226 

2120 

0.847 

3.269 

21.2 

24.9 

2 

1.56 

2960 

17.8 

234 

2141 

0.848 

3.258 

19.4 

23.0 

3 

1.64 

2945 

17.0 

225 

2110 

0.846 

3.213 

14.2 

19.6 

4 

1.42 

2978 

17.4 

233 

2196 

0.852 

3.219 

14.6 

18.9 

5 

1.50 

2963 

17.0 

207 

2164 

0.850 

D.209 

16.2 

19.8 

6 

1.66 

2940 

17.1 

229 

2104 

0.846 

3.130 

7.01 

11.45 

7 

1.72 

2935 

17.5 

220 

2082 

0.844 

3.128 

9.86 

11.65 

8 

1.41 

2978 

17.1 

209 

2199 

0.852 

3.262 

18.8 

24.6 

9 

1.51 

2963 

16.9 

214 

2159 

0.849 

3.262 

20.7 

24.8 

10 

1.75 

2930 

17.7 

225 

2072 

0.844 

0.290 

20.4 

27.7 

11 

1.67 

2942 

17.7 

224 

2100 

0.845 

0.296 

21.7 

28.0 

ethanol : 

12  1.18 

11 

17.8 

223 

n 

II 

i 

0.809 

14.2 

50.0 

13 

1.10 

11 

17.4 

222 

n 

II 

I 

0.869 

16.2 

54.4 

14 

1.17 

11 

17.8 

225 

it 

II 

1 

0.834 

17.4 

51.4 

15 

1.12 

vv 

17.4 

211 

ii 

II 

1 

0.859 

19.4 

55.2 

17 

1.25 

II 

18.0 

209 

ti 

n 

0.872 

17.7 

56.1 

18 

1.17 

n 

17.7 

228 

ii 

N 

0.869 

18.0 

53.4 

19 

1.13 

ri 

17.6 

224 

ii 

n 

0.716 

16.2 

43.8 

20 

1.25 

»i 

18.0 

224 

n 

it 

0.713 

16.2 

43.5 

22 

1.41 

2978 

18.1 

236 

2199 

0.852 

0.570 

14.6 

33.1 

23 

1.37 

11 

17.9 

225 

II 

II 

0.598 

15.8 

35.9 

25 

1.48 

II 

18.1 

230 

II 

11 

0.535 

11.0 

31.4 

26 

1.31 

II 

18.2 

237 

II 

II 

0.525 

10.4 

30.2 

28 

1.54 

II 

18.0 

238 

N 

II 

0.341 

10.3 

19.0 

29 

1.53 

11 

17.8 

225 

If 

n 

0.369 

9.55 

21.4 

30 

1.59 

n 

18.2 

231 

11 

H 

0.318 

10.2 

17.9 
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A  comparison  of  the  calculations  with  the  experimental  data  for  water 
coolant  is  given  in  Figure  8.  The  calculated  liquid  film  lengths  are  an  average  of 
25%  longer  than  measured.  Since  convection  accounted  for  62%  of  the  total 
evaporation  rate,  multiplying  the  convective  heat  flux  by  a  factor  of  1.4  would  bring 
the  calculations  into  best  agreement  with  the  measurements.  This  corresponds  to  a 
free  stream  turbulence  of  9%,  which  is  a  reasonable  assumption. 

In  all  of  the  ethanol  tests  the  liquid  coolant  flow  rate  (T)  greatly  exceeded 
Knuth’s  critical  value  for  the  formation  of  large  waves  (T^),  given  in  the  coolant 
table  above.  Such  large  waves  would  explain  why  the  measured  liquid  film  lengths 
were  an  average  of  2.75  times  shorter  than  calculated,  which  compares  favorably 
with  Kinney’s  results  for  the  mass  loss  rate  with  large  waves  present.  For  simplicity, 
die  gas  properties  at  an  O/F  of  1.41  were  used  in  all  of  the  ethanol  calculations. 

Katto's  correlation  for  radiant  burnout  is  tested  in  the  program  and  predicted 
liquid  film  lengths  10%  shorter  than  listed  above  with  water  coolant.  However,  with 
ethanol  coolant  the  burnout  condition  predicted  liquid  film  lengths  of  less  than 
1  cm,  which  is  unreasonable.  Therefore  the  burnout  prediction  was  not  deemed 
trustworthy. 
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4.1.4  Warner  and  Emmons  H 2/air  rocket  (Reference  34) 

These  tests  were  made  in  a  79  8  mm  diameter  rocket  chamber  burning  air 
and  H2.  The  liquid  film  cooling  test  section  was  276  mm  downstream  of  the 
injector.  The  data  given  in  their  Figure  5  was  analyzed.  The  gas  was  assumed  due 
to  the  stoichiometric  reaction  of  air  and  Ha,  giving  0.653  Na  and  0.347  HaO  molar 
fractions,  and  a  molecular  weight  of  24.5.  A  temperature  of  2222  K  and  a  pressure 
of  34  atm  were  reported. 

Since  the  mass  flow  of  gas  (Gcj,)  was  not  given,  it  was  back-calculated  from 
the  stated  Reynold’s  number  of  68,000.  At  the  free  stream  temperature  die  species 
viscosities  were  6.626  10’ 6  kg/m-s  for  N2  and  7.186  10'6  kg/m-s  for  HaO,  Using 
Wilke’s  method,  the  free  stream  gas  viscosity  was  6.806  10‘5  kg/m-s.  Using  this 
value:  Gcjj  -  58.03  kg/s-m2. 

The  gas  properties,  evaluated  at  a  mean  temperature  of  1300  K,  are  listed 

below: 

Cp  Pr  v 

(J/kq-K) _ fl0~6kq/m-sl_ 

Na  1226  0.705  4.6B0 

Ha0  2494  0.90  4.5B4 

mixture  1549  0.773  4.670 

The  property  values  were  determined  in  the  same  manner  as  for  Morrell’s  data. 


Three  different  coolants  were  tested.  The  properties  of  each  at  the  chamber 
pressure  are  given  below: 


coolant 

Tv 

(K) 

i  Cgi  nv 

f  10— a 

i  kail ka-K i  1 

kg_l 

m-sl 

fiO-’Hl 

1  mL 

P 1  Pv 
f  1 

rcr 

'kg_l 

s-mi 

1  m3  1 

water 

514 

1.76  4779  17.1 

111.0 

2B.4 

813  17.0 

0.266 

ethanol 

480 

0.440  5126  14.0 

83.0 

4.18 

564  50.3 

0.239 

ammonia 

345 

0.937  4103  11.9 

79.4 

10.25 

523  25.2 

( 

3.180 

The  coolant  properties  were  determined  as  In  Morrell’s  tests,  with  the  exceptions 
that  X,  Cpj,  p\,  and  pv  for  ammonia  were  scaled  from  the  values  at  325  K  in  Table 
3-209  of  Perry  (Reference  56). 


With  these  property  values,  the  liquid  film  length  and  average  evaporation 
rates  for  each  test  were  calculated  at  a  coolant  flow  (T)  of  0.08  kg/s-m,  giving: 
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Coolant 

Lc 

fern) 

®vcalc 
fkcr/s-m* 1 

mvexp 
fkcr/s-m2 1 

water 

29.05 

0.275 

0.439 

NH3 

14.25 

0.561 

1.159 

ethanol 

16.85 

0.475 

0.886 

Both  the  ammonia  and  ethanol  data  are  somewhat  questionable  since  they 
extrapolate  to  a  positive  coolant  flow  at  zero  film  cooled  Length,  indicating  an  initial 
mass  loss  at  Injection.  The  ammonia  data  also  showed  significant  scatter.  The  flow 
rates  in  the  ammonia  tests  were  centered  about  Knuth’s  critical  value,  which  may 
explain  the  scatter. 

The  calculated  evaporation  rate  for  water  coolant  is  37%  lower  than 
measured.  Radiation  accounted  for  62.8%  of  the  calculated  evaporation  rate, 
making  the  calculations  very  sensitive  to  the  assumed  gas  temperature.  Given  the 
uncertainties  in  die  test  conditions,  it  is  futile  to  speculate  at  length  as  to  the 
discrepancy. 
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4.1.5  R.C.  Kesselring,  et  al  OFa/BaHa  rocket  test  (Reference 53) 


These  tests  were  performed  in  a  rocket  engine  using  diborane  fuel,  which  is 
uncommon  and  very  expensive.  The  combustion  chamber  was  a  thin  niclde  shell  of 
9.45  cm  diameter,  used  for  short  duration  tests.  The  free  stream  gas  conditions 
were  fairly  constant  for  all  runs,  so  that  the  average  properties  below  were  used  in 
analyzing  all  runs: 

Tg=  3900  K,P= 6.8  atm,  Gch= 150.7  kg/s. m3  . 

The  NASA/Lewis  Chemical  Equilibrium  code  gave  the  following  free  stream  mole 
fractions  and  properties: 

0.537  HF,  0.225  BOF,  0.063  Ha  O,  0.036  Ha ,  0.031  OH,  Mg=25-6 
and  at  Tmean  =  2200  K: 

Cpg=  3665  J/kg-K,  ng =8.37  10'*  kg/m-s,  Pr=0.80 

The  properties  of  saturated  BaHe  coolant  at  6.94  atm  from  Kit  and  Evered 
(Reference  58)  are: 

Tv=  243.8  K,  Mc=  27.69,  X  =  3-72  10s  J/kg,  Cpl=  3235  J/kg-K 

(i\=  3  6  lO^kg/ims,  p\  =  341  kg/ms ,  pv-  29.4kg/m*,a  =  7.1610"3N/m 

Kit  and  Evered  give  X  -  5.173  10*  J/kg  at  Tjj=  180.04  K  at  1  atm,  which  compares 
well  with  Perry’s  handbook  (Reference  56).  Scaling  X  as  suggested  in  Reid 
(Reference  55)  to  243.8  K,  using  Tcrit=  289.7  K,  gives  the  value  above.  The  vapor 
viscosity  was  not  found.  However,  due  to  the  very  low  liquid  viscosity,  the  critical 
flow  rate  for  forming  large  waves,  rcp  is  estimated  to  exceed  100  kg/m.s,  making 
such  waves  unlikely. 

Thermocouples  were  located  on  the  outside  of  the  thin  chamber.  A  steady 
temperature  rise  during  a  run  was  indicative  of  the  absence  of  liquid  film  on  the 
wall  at  that  axial  position.  Where  the  liquid  film  did  exist  the  wall  temperature 
actually  decreased  during  a  firing,  due  to  the  low  saturation  temperature  of 
diborane.  Unfortunately,  the  thermocouples  were  spaced  2  inches  apart,  making  it 
difficult  to  determine  the  liquid  film  lengths  accurately.  The  best  approximations 
are  given  below: 
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TEST 

T  fka/m-sl 

Lc  fin) 

1 

0.583 

3<  Lc  <5 

3 

0.445 

1<  LC  <3 

4 

0.451 

1<  LC  <3 

€ 

0.191 

1<  Lc  <3 

9 

0.306 

1<  Lc  <3 

10 

0.532  . 

“  3 

11 

0.435 

1<  LC  <3 

Surprisingly,  die  authors  assume  that  the  liquid  is  immediately  evaporated, 
based  upon  calculations  of  the  normally  expected  heat  flux  without  transpiration. 
They  interpret  the  thermocouple  data  as  meaning  that  this  vapor  refrains  from 
mixing  with  the  hot  gases  until  some  arbitrary  distance  downstream. 

No  information  on  the  emittance  of  the  combustion  products  was  available, 
so  the  emittance  was  treated  as  an  adjustable  parameter.  Assuming  8%  free  stream 
turbulence,  the  best  comparison  between  the  data  and  calculations  was  obtained 
with  an  assumed  emittance  of  17%,  as  shown  in  Figure  9.  Radiation  accounted  for 
90%  of  die  heat  flux  in  this  case.  The  calculated  liquid  film  length  is  somewhat 
insensitive  to  the  assumed  emittance.  For  example,  halving  the  emittance  to  8.5% 
increases  the  calculated  liquid  film  length  by  58%,  whereas  an  82%  increase  might 
be  expected.  The  reason  is  that  the  vapor  liberated  by  radiation  greatly  decreases 
die  convective  heat  flux  through  transpiration.  If  the  radiative  flux  decreases, 
convection  increases,  keeping  the  evaporation  rate  somewhat  constant. 

Given  the  numerous  species  in  the  products,  one  might  expect  a  larger 
emittance,  however  the  optical  density  of  64  atm-cm  was  fairly  low.  To  critically 
judge  the  analysis  for  this  test  requires  an  independent  determination  of  the  gas 
emittance. 
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42  GASEOUS  FILM  COOLING 


4.2.1  Ewen  and  Rousar  (References  49-51) 


The  goal  in  all  three  test  programs  was  to  study  gaseous  film  cooling  in  a 
converging-diverging  nozzle,  for  application  to  rocket  engines.  The  authors 
analyzed  their  data  in  terms  of  an  empirical  model.  Although  their  model  is  similar 
to  the  differential  entrainment  formulation  of  Section  3.4.1,  it  is  not  based  upon 
standard  correlations.  For  this  reason  only  their  raw  data  was  used. 


In  the  first  two  reports  (1972  and  1973)  testing  was  done  on  a  flow  bench  in 
a  laboratory,  using  heated  nitrogen  for  the  free  stream  flow  and  either  hydrogen  or 
nitrogen  coolant.  The  free  stream  turbulence  was  measured  to  be  4.1%  for  all  tests. 
The  free  stream  flow  was  folly-developed  upon  contacting  the  coolant  injector, 
leading  to  a  difficulty  in  specifying  the  injection  position  from  the  "leading  edge",  a 
problem  discussed  later.  Only  those  tests  with  a  coolant  injection  velocity  less  than 
die  free  stream  velocity  were  analyzed  (to  avoid  jet  effects),  although  all  of  the  data 
was  qualitatively  similar.  The  test  conditions  are  given  below: 


Test 

chamber 

Gch 

fka/s»m2 1 

coolant 

MC 

fkq/s»ml 

Tc 

m 

Aet 
_ (41 

1972: 

3 

30  * -1R 

630 

763 

H2 

0.189 

294 

38 

5A 

M  II 

629 

764 

Ha 

0.185 

285 

28 

8 

II  II 

636 

756 

Ha 

0.122 

111 

6 

11A 

conical 

631 

761 

h2 

0.185 

292 

0 

1973: 

101L 

30*-lR 

581 

820 

n2 

0.355 

341 

0 

101H 

ii  n 

577 

822 

tJ2 

0.528 

320 

0 

102 

30’-2R 

528 

880 

H2 

0.170 

270 

16 

103A 

15--2R 

591 

816 

0.173 

275 

18 

104AL 

ii  ii 

599 

815 

n2 

0.360 

345 

0 

104AH 

ii  ii 

596 

817 

N2 

0.677 

327 

0 

105A 

con .  mod 

597 

817 

h2 

0.157 

270 

0 

For  test  104AL,  the  free  stream  temperature  was  stated  as  915’  K  in  Rousar 
and  E wen's  report,  apparently  as  a  misprint.  Also,  for  tests  102,  103A,  and  105A, 
the  stated  coolant  injection  temperatures,  Tc,  were  adjusted  to  be  consistent  with 
the  axial  thermocouple  readings. 


In  the  last  report  (1977),  measurements  were  made  in  an  02(g)/H2  rocket 
engine  using  hydrogen,  nitrogen,  and  helium  coolants.  All  tests  used  the  same 
chamber  (P/N  1182134).  The  coolant  injector  was  178  mm  downstream  of  the 


43 


AEDC-TR-91-1 


propellant  injector.  The  test  conditions  are  given  below: 


Test 

O/F 

Gch 

(ka/s.m2 1 

Si 

coolant 

Mc 

(ka/s.m) 

Tc 

(K) 

Aet 

m 

102 

3-89 

242.6 

2670 

H2 

0.587 

290 

0. 

103 

3-81 

228.5 

2640 

H2 

0.504 

290 

0. 

105 

5-75 

241.3 

2990 

H2 

0.478 

290 

0. 

110 

7.71 

263.1 

3080 

H2 

0.414 

290 

0. 

111 

7.99 

260.6 

3080 

H2 

0.362 

290 

0. 

114 

6.27 

237.5 

3030 

H2 

0.345 

290 

0. 

115 

5.93 

234.9 

3005 

H2 

0.408 

290 

0. 

116 

7.60 

246.4 

3075 

H2 

0.309 

290 

0. 

117 

4.00 

219.5 

2720 

H2 

0.391 

290 

0. 

119 

7.76 

283.7 

3080 

He 

0.515 

310 

0. 

120 

7.46 

270.8 

3075 

He 

0.469 

310 

0. 

122 

7.75 

299.1 

3080 

He 

0.595 

310  ' 

0. 

The  coolant  flow  rates  (W^  were  not  given  directly  and  were  back-calculated 
from  the  stated  total  mass  flow  (Wp),  die  coolant  to  fuel  mass  flow  ratio  (Wc/Wp), 
and  the  O/F  ratios  as: 


Wc  =  (WC/WF)  Wt/(1+0/F) 

The  Mc  values  tabulated  above  are  die  coolant  mass  flows  (W^  per  circumference. 

In  all  of  the  tests  die  products  contained  excess  hydrogen.  From  the 
reaction  stoichiometry,  the  mass  fraction  of  HzO  in  the  products  is  given  by: 

®h20  =  9/8  (0/P)/ [1+ (0/P)] 

the  remainder  being  H3.  The  flame  temperatures,  Tg,  listed  are  from  a  plot  of 
somewhat  conflicting  values  from  a  number  of  literature  sources,  mostly  at  1  atm 
pressure. 

The  specific  heats  of  the  various  gas  species  were  calculated  as  (Reference  57): 

Cp  |  N  =  1415  -  2.88  10®  K/T  +  53  46  10«  K2  /T2  J/kg.  K 
CpiH2  =  12,047  +  2.176K-1  T  +  3U80KK/T*  J/kg.K 
Cp|  H*o  =  4615  - 1-03  10B ¥?• /T*4  +  9.68  10SK/T  J/kg-K 
Cp| He  =  5238  J/kg.K 

These  specific  heats  were  weighted  by  the  species  mass  fractions  to 
determine  the  product  specific  heats. 
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The  viscosity  of  each  species  was  calculated  from  the  Chapman-Enskog 
constants.  The  mixture  weighting  method  of  Wilke  was  used  to  determine  the 
product 'viscosity  (Reference  55): 

/I  In  -  1.23  10T«  JT  kg/m.5 
*jjH2  =6.711tr7JT  kg/m.s 
^Ih'o  =  822  HT'JT  kg/m.s 

The  parameters  describing  the  chamber  contours,  discussed  in  Section  3.6, 
are  tabulated  below: 


chamber 

Dch 

(mm) 

ll 

(mini 

rl 

fmm) 

*C 

(deal 

(mm) 

r2 

(mm) 

*d 

(deal 

30 ’ -1R 

31. 

88.9 

15.5 

30. 

15.5 

7.8 

15. 

30,-2R 

31. 

88.9 

31. 

30. 

15.5 

15.5 

15. 

15  * -2R 

31. 

88.9 

31. 

15. 

15.5 

15.5 

15. 

conical 

31. 

7.9 

12.7 

4.4 

15.5 

7.8 

15. 

(modified) 

31. 

0. 

12.7 

4.4 

15.5 

7.8 

15. 

1182134 

99.6 

106.7 

51.8 

30. 

52.1 

28.7 

15. 

The  results  of  each  calculation  are  compared  with  the  data  in  Figure  10.  The 
simplest  geometry  was  test  11A,  with  a  gradually  converging  conical  chamber.  For 
comparison,  the  calculations  are  shown  both  with  the  differential  entrainment 
formulation  (Section  3.4.1)  and  with  the  simpler  integral  correlation  (Section  3.0) 
(based  upon  the  starting  conditions).  The  injection  point  was  assumed  at  xjn=0  for 
both.  The  differential  formulation  gives  the  best  comparison  up  to  the  throat. 
Downstream  of  the  throat  the  wall  temperature  decreases  much  fester  than 
predicted  by  the  model,  a  difficulty  with  all  of  the  calculations  which  is  discussed 
later. 


The  comparison  is  not  as  good  for  the  cylindrical  chamber  data,  such  as  for 
test  3.  The  wall  temperature  increases  very  rapidly  downstream  of  the  converging 
turn  (at  88.9  mm).  The  calculated  value  might  be  expected  to  rise  due  to  die 
increase  in  the  free  stream  mass  flow  as  the  chamber  area  decreases.  However,  as 
shown  by  the  solid  curve,  this  is  not  the  case.  The  reason  is  that  the  contraction 
also  increases  the  boundary  layer  mass  flow  per  circumference,  as  discussed  in 
Section  3.4.2,  and  the  two  effects  tend  to  cancel.  Disabling  this  circumferential 
correction  term  yields  a  better  comparison,  as  shown  by  the  dashed  curve. 
However,  such  fortuitous  comparison  is  not  a  suitable  reason  for  discounting  a 
physical  effect,  so  the  circumferential  term  is  left  in  the  model. 
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As  in  most  modeling  efforts  the  comparison  can  be  improved  by  judicious 
selection  of  input  parameters,  as  shown  in  the  second  figure.  To  force  comparison 
in  the  straight  chamber  section,  a  downstream  injection  point  (x^)  of  20  mm  was 
used.  'While  arbitrary,  this  could  be  defended  by  arguing  that  in  fully-developed 
flow  the  injection  point  is  considered  an  infinite  distance  from  the  "leading  edge". 
However,  with  the  exception  of  tests  3  and  5A,  xjn=0  was  always  the  optimum 
choice.  This  is  supported  by  heat  transfer  data  taken  by  Ewen  and  Rousar  which 
show  the  initially  high  convective  heat  transfer  rates  of  a  leading  edge  at  the 
injection  point  Apparently  the  coolant  injection  disrupts  the  flow  such  that  a  new 
boundary  layer  grows  from  the  injection  point. 

To  force  comparison  downstream  of  the  converging  turn  it  was  necessary  to 
assume  an  increase  in  the  free  stream  turbulence  due  to  the  turn,  Ae*,  of  38%  for 
test  3-  Similar  correction  factors  are  necessary  with  all  turbulence  models 
(Reference  59).  The  values  required  for  each  test  are  tabulated  above. 


Such  a  correction  is  defensible  only  if  it  is  not  arbitrary.  This  extra  "turning 
turbulence”  can  be  correlated  by  a  "centrifugal  parameter": 


Fc=  ip  =  Si 

c  g  c  rt  rlPg 


where  rt  is  the  radius  of  the  converging  turn,  pQ  and  p g  are  the  coolant  and  gas 
densities  at  the  turn.  Note  that  Fc  has  dimensions.  For  all  of  the  1977  tests  in  the 
rocket  engine,  the  correction  was  negligible.  Given  the  low  densities  in  most  rocket 
chambers,  this  is  generally  true. 


As  shown  in  Figure  11,  the  correlation  between  Aet  and  Fc  can  be  fit  by  a 
function: 


Aet  =  5.26  ( [ 1+ (0 . 029  Pc)®]1/*  -  1)  for  Pg  >  PC 
=  o.  for  pg  <  pc 

where:  Fc  <>  kN/ms 

The  aberrant  point  is  due  to  test  8  of  1972,  for  which  it  was  impossible  to  obtain 
any  satisfactory  comparisons. 

It  was  not  possible  to  obtain  a  reasonable  comparison  with  the  data  of  tests 
119-122  which  used  helium  coolant.  The  rate  of  temperature  rise  with  non-reactive 
helium  coolant  was  overpredicted  by  a  factor  of  1.7,  raising  the  question  that 
combustion  effects  are  important  and  that  the  favorable  comparisons  with  hydrogen 
coolant  in  the  rocket  were  merely  fortuitous. 
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For  instance,  the  free  stream  temperature  dose  to  the  wall  might  have  been 
considerably  less  than  that  used  in  the  calculations  above.  Since  many  injectors  provide 
a  fuel  rich  zone  near  the  wall,  this  is  plausible.  Such  an  assumption  would  make  the  He 
coolant  calculations  compare  more  favorably  with  the  data.  Normally  this  assumption 
would  also  lead  to  a  significant  under-prediction  of  the  wall  temperatures  with  H2  coolant. 
However,  considering  combustion  of  the  H2  coolant  in  the  boundary  layer  could  restore 
the  wall  temperatures  to  the  measured  values.  Data  with  other  coolants  is  needed  to  ascer¬ 
tain  whether  combustion  effects  could  explain  the  discrepancy. 

In  all  of  the  calculations  it  was  not  possible  to  account  for  die  rapid  decrease 
of  wall  temperature  downstream  of  the  throat.  Since  the  wall  temperature 
decreased  even  faster  than  the  free  stream  gas  static  temperature,  increased  mixing 
with  the  free  stream  gases  could  not  account  for  this  rapid  decrease.  Instead,  It 
might  be  attributed  to  the  static  temperature  change  of  the  boundary  layer  gases  as 
they  accelerate.  Since  die  present  formulation  lumps  together  all  of  the  boundary 
layer  gases,  it  would  be  difficult  to  incorporate  such  an  effect  into  the  model.  A 
turbulent  boundary  layer  model,  which  considers  the  velocity  profile  across  the 
boundary  layer,  would  be  better  suited  to  this  task. 

All  of  the  calculated  results  assumed  a  recovery  factor  of  zero,  meaning  that 
the  entrained  free  stream  gas  is  at  the  static  temperature.  However,  upstream  of  the 
nozzle  throat  there  is  little  difference  between  the  static  and  stagnation 
temperatures,  so  that  the  choice  of  recovery  factor  is  not  significant.  Downstream 
of  the  throat  the  present  model  is  unable  to  accurately  predict  the  wall  temperature, 
regardless  of  the  choice  of  recovery  factor. 
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4 22  G.R.  Kinney,  et  al  (Reference  17) 

Hie  data  in  this  report  was  reviewed  in  Section  4.1.1  pertaining  to  calculation 
of  Liquid  film  lengths,  which  was  the  focus  of  their  study.  The  authors  do,  however, 
give  one  plot  of  temperatures  downstream  of  the  dryout  point  in  their  Figure  4. 
This  data  is  useful  in  determining  how  to  apply  foe  gaseous  film  cooling  analysis  to 
foe  special  case  of  liquid  film  cooling.  The  property  values  are  from  Section  4.1.1 . 

The  start  of  foe  liquid  film  was  40  inches  from  the  tube  entrance,  and  foe 
liquid  film  was  20  inches  long.  As  argued  in  the  previous  analysis  of  Ewen  and 
Rousar’s  data,  the  gaseous  coolant  injection  tends  to  disturb  the  boundary  layer 
such  that  foe  injection  point  acts  effectively  as  a  "leading  edge".  The  simplest 
assumption  is  that  all  of  foe  vapor  is  injected  at  foe  termination  of  foe  liquid  film 
(xjn=0).  The  calculated  wall  temperatures  compare  favorably  with  the  data,  as 
shown  in  Figure  12.  This  is  compared  with  an  alternate  assumption  that  the 
effective  leading  edge  is  in  foe  center  of  foe  liquid  film  (xjn  =  254  mm),  which  gives 
very  poor  comparison  with  the  data.  Therefore,  as  was  found  with  Ewen  and 
Rousar’s  data,  taking  Xj„  =  0  generally  gives  the  best  comparisons  with  existing 
data  for  all  cases. 
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5.0  CONCLUSIONS 

A*  simple  one-dimensional  model  gives  satisfactory  comparison  with  existing 
data  for  liquid  film  lengths  in  rocket  engines.  Convective  and  radiant  heat  transfer 
must  bodi  be  considered.  At  high  coolant  flow  rates  large  waves  on  the  film  can 
decrease  the  coolant  efficiency.  Knuth’s  correlation  Is  useful  in  predicting  this 
transition.  Radiant  burnout  of  the  liquid  film  is  possible,  however  the  use  of 
existing  burnout  correlations  is  questionable. 

Downstream  of  die  liquid  film  a  standard  gaseous  film  cooling  correlation, 
modified  to  a  differential  form,  satisfactorily  predicts  the  wall  temperatures 
upstream  of  the  throat  as  the  vapor  mixes  with  the  free  stream  gas.  A  correction 
factor,  correlated  in  terms  of  a  "centrifugal  parameter",  is  needed  to  account  for 
increased  mixing  due  to  the  converging  turn.  However,  in  most  rocket  engines  it  is 
not  significant  Downstream  of  the  throat  the  wall  temperature  drops  very  rapidly 
due  to  acceleration  of  the  boundary  layer  gases,  which  is  not  predicted  by  the 
present  model. 
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NOMENCLATURE 


A  =  local  cross-sectional  area  of  chamber 
A w  =  absorptivity  of  chamber  walls 
Cf  -  skin  friction  factor  “  2rw/pgUg2 
Cp  =  specific  heat  per  mass 

D  =  diameter  of  combustion  chamber  at  position  x 

et  =  free  stream  turbulence  intensity  fraction 

Aet=  increase  in  turbulence  due  to  converging  turn 

F  =  blowing  parameter  =  nty/G 

Gch  =  free  stream  gas  mass  flow  per  area  =  pUg 

Gmean  =  mass  A°w  at  Tmean  relative  to  liquid  surface  = 

h  =  convective  heat  transfer  coeff. 

H  =  F/St 

K  =  G^g0,26/  Me1 ,26 
Kg  =  thermal  conductivity  of  free  stream  gas 
“  correction  factor  for  molecular  weight 
Kt  =  correction  factor  for  turbulence 
L  =  optical  path  length 
Lg^-  =  average  optical  pad)  length 
mv  =  total  liquid  evaporation  rate  per  surface  area 
M  =  molecular  weight  or 

Mjj[  =  boundary  layer  mass  flow  rate  per  circumference 

Mc  =  gaseous  coolant  mass  flow  rate  per  circumference 

n  =  no.  of  moles 

Nc  =  mole  fraction  of  CO2 

Nw  =  mole  fraction  of  water  vapor 

Fr  =  Prandtl  no  of  gas  e  PgCpg/Kg 

P  =  absolute  pressure 

Q  =  heat  flow  per  area 

r  =  radius  of  convergence  arc  in  nozzle 

ReD  =  Reynold's  no  based  on  diameter  =  GD//jg 

St  =  Stanton  no  m  h/(G  Cpg) 

T  =  absolute  temperature 
Tc  =  temperature  of  liquid  or  gaseous  coolant 
Tv  =  saturation  temperature  of  coolant 
hT  -  Tg-Ty 
U  -  axial  velocity 


PmeanCUg’Ul) 
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x  =  distance  from  injection  point  (along  contour) 

Xg  =  x  corrected  for  developing  pipe  flow 
X  =  dimensionless  distance  =  K  x 
y  =  distance  from  wall 

Greek  Symbols:, 

7  =  Cp/Cy  =  specific  heat  ratio 

S  =  boundary  layer  thickness 

e  =  gas  emissivity 

j)  -  film  cooling  effectiveness 

T  -  liquid  coolant  mass  flow  rate  per  circumference 

X  =  latent  heat  of  vaporization  of  coolant 

X*  =  X  +  CpjCIVTc) 

fi  =  dynamic  viscosity 

p  =  mass  density 

Popt  =  optical  density  =  P  L 

o  =  surface  tension  of  coolant  or  Stephan’s  const. 

subscripts: 

aw  -  adiabatic  wall 

bl  -  boundary  layer 

bo  -  burnout  point 

c  -  coolant,  (liquid  or  gaseous) 

ch  -  values  in  cylindrical  chamber  section 

conv  -  convective 

cr  -  at  transition  to  "large  waves" 

g  -  free  stream  gas 

1  •  liquid  coolant 

mean  -  mean  temperature  between  free  stream  gas  and  liquid 
o  -  for  "dry-wall"  conditions,  without  transpiration 
rad  -  due  to  radiation 
up  -  upstream 
v  -  vapor 

w  -  evaluated  at  wall 
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